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ABSTRACT 

We present derived stellar and disc parameters for a sample of Taurus brown 
dwarfs both with and without evidence of an associated disc. These parameters have 
been derived using an online fitting tool (http://bd-server.astro.ex.ac.uk/), 
which includes a statistically robust derivation of uncertainties, an indication of pa- 
rameter degeneracies, and a complete treatment of the input photometric and spec- 
troscopic observations. 

The observations of the Taurus members with indications of disc presence have 
been fitted using a grid of theoretical models including detailed treatments of physi- 
cal processes accepted for higher mass stars, such as dust sublimation, and a simple 
treatment of the accretion flux. This grid of models has been designed to test the va- 
lidity of the adopted physical mechanisms, but we have also constructed models using 
parameterisation, for example semi-empirical dust sublimation radii, for users solely 
interested in parameter derivation and the quality of the fit. 

The parameters derived for the naked and disc brown dwarf systems are largely 
consistent with literature observations. However, our inner disc edge locations are 
consistently closer to the star than previous results and we also derive elevated ac- 
cretion rates over non-SED based accretion rate derivations. For inner edge locations 
we attribute these differences to the detailed modelling we have performed of the disc 
structure, particularly at the crucial inner edge where departures in geometry from 
the often adopted vertical wall due to dust sublimation (and therefore accretion flux) 
can compensate for temperature (and therefore distance) changes to the inner edge 
of the dust disc. In the case of the elevated derived accretion rates, in some cases, 
this may be caused by the intrinsic stellar luminosities of the targets exceeding that 
predicted by the isochrones we have adopted. 

Key words: starsievolution - stars:formation - stars: pre-main-sequence - tech- 
niques: photometric - catalogues - (stars) Hertzsprung-Russell H-R diagram 



1 INTRODUCTION 

In studies of star formation and evolution spectral energy 
distribution (SED) fitting is often used to derive parameters 
for stars and their circumstellar discs. Fitting tools, such as 
[Robitaille et al.] ( |2006[ ), and the increasing availability of 
observational data have led to a rapid increase in the num- 
ber of studies using SED fitting. As observations of fainter 
lower-mass stars (LMS) and brown dwarfs (BDs) have be- 



2003 



come feasible (see for example Jayawardhana et al, 

Luhman et "ari[2005l [20081 |Guieu et al.||2007| |Monin et al.| 
2010 Olofsson et al. 2010) such SED fitting techniques have 
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been extended into the lower mass regime. These studies 
have generally been limited to small numbers of members of 
the Taurus star forming region (SFR) due to its proximity. 

For higher mass stars extensive studies have been per- 
formed using detailed treatments of the important physi- 
cal processes, for example radiative transfer, dust sublima- 
tion and vertical hydrostatic equilibrium (see for example 
Tannirkulam et al.|[2008i Walker et al.||2004[ ). Such models 



are computationally expensive and as such are usually used 
to model individual systems. The online fitting tool of |Ro-| 



bitaille et al. ( 2006 1 has allowed much greater access to such 
models via pre-computed 'grids' of radiative transfer simu- 
lations covering a selected parameter space, for example a 
range of masses and ages. However, the young-stellar objects 



2 N.J. Mayne et al. 



(YSOs) modelled by |Robitaille et al.| p006| ) do not extend 
into the BD mass regime. 

In the low mass regime, studies of populations of ob- 
jects are generally performed using parameterisations of 
physical mechanisms rendering the problem computation- 
ally tractable, x^ fitting is used to derive parameters from 
a grid of model SEDs computed over ranges of parameters 



selected by the user (e.g. Bouy et al. 2008 Monin et al, 
2010 Qlofsson et al.pOlO l. For example, the inner bound- 
ary of the associated dust disc can be placed at a range of 
positions, the resulting SEDs derived, and the best fitting 
model selected, leading to an estimate of the inner edge ra- 
dius. However, it is known that the fiux from the inner edge 
is a complicated function of its geometry, composition and 
temperature, all of which can change as a function of the 
stellar properties (see discussion in [Mayne fc Harries|20l6| ). 
A clear correlation between the observed inner edge location 
and the dust sublimation radius has been found for higher 
mass stars ( [Monnier fc Millan-G abct"2002). Our approach 
in this work is to constrain our modelling of BD and disc 
(BDD, hereafter) systems by adopting the physical mecha- 
nisms that have proved successful on their higher mass coun- 
terparts such as Classical T Tauri stars and Herbig Ae stars: 
(i) we constrain our disc structures so that they are in ver- 
tical hydrostatic equilibrium, (ii) The inner edge of the gas 
disc is truncated at the co-rotation radius (as such it is a 



function of stellar mass and rotational period, Camenzind 



19901 |Koenigl|[T99Tl [MuzeroUe et~aL|[2003l [Mohanty fc^iiu 



2008 1. (iii) The location of the dust disc inner radius is de- 
fined dust sublimation ( Tannirkulam et al.||2008 l. (iv) The 
photosphere luminosity is enhanced by accretion luminosity 



from material impacting on the BD surface (Bouvier et al. 
1995 Herbst et al.|2007 l. Note that although we are appar- 
ently increasing the physical complexity of the model, we are 
in fact reducing the number of free parameters. As a conse- 
quence of this adequate fits to the SEDs act to corroborate 
the physical assumptions that underly our models, whilst 
inadequate fits lead us to question them. In contrast allow- 
ing (for example) the the disc structure to vary arbitrarily 
will necessarily result in an improved spectral fit, but per- 
haps at the cost of physical insight. Applying such models 
to the studies of populations then requires sufficient com- 
puting resource, or an available repository of pre-computed 
models. 

Comparison of data and models is further complicated 
by temporal variability. Simulations of the dynamics of cir- 
cumstellar discs have suggested that several mechanisms 
can lead to structural variations over a range of timescales, 
some as short as a few days (|Min et al.||2009{ |Acreman| 



et al. 20101. Combined with variations in emissions from 



the stellar photosphere from processes such as accretion 
( [Bouvier et al.||1995| [Herbst et "aL]|2007| l , this leads to sig- 
nificant variations in the observable properties. However, 
non-contemporaneous photometric data sets are often com- 
bined and fitted to derived SEDs, without consideration 
of such temporal variability. Additionally, previous studies 
have found that the optical photometry of BDs and LMS, 
which is vital for constraining the stellar parameters, are 
poorly represented by simulated stellar atmospheres. Fur- 
thermore, as the slope of the SED across the optical regime 
is much steeper than for higher mass systems, these measure- 
ments dominate the y^ value for a given fit and so are often 



omitted from the fit (see for example Guieu et al. 2007 1. Dis- 
crepancies in the optical photometry, however, are often the 
result of converting magnitudes to a standard photometric 
system (to match the models), or indeed converting mag- 
nitudes to monochromatic fiuxes, using simple transforma- 
tions or calibrations. Such transformations are usually only 
accurate for main-sequence (MS) stars, and their use will 
introduce additional uncertainty in the observed BD fluxes. 

We have used the TORUS radiative transfer code to sim- 
ulate a large number of BDD systems covering a range of 
input parameters and adopting mechanisms accepted for 
higher mass systems. These simulations include a represen- 
tation of dust sublimation, vertical hydrostatic equilibrium, 
accretion flux and radiative transfer. These models have 
been incorporated into an online fitting tool freely accessible 
to the community at http://bd-server.astro.ex.ac.uk/. 
This fitting tool includes an accurate treatment of the obser- 
vational photometry so that it can be compared consistently 
with the simulated observations. As a proof of concept, and 
to aid development of this tool, we have fitted a sample of 
Taurus BD and BDD observations collected from the litera- 
ture, the results of which we present in this paper. We have 
explored differences in our derived parameters from those in 
the literature and discuss some specific examples. 

This paper is structured as follows. Section [2| intro- 
duces the adopted data-set for the sample of Taurus BDs 
and BDDs. The detail of the transformations and genera- 
tion of uncertainties (the latter to be applied to the model 
during fitting) are mentioned in this section and explained in 
full in Appendix [A] The models used, including the physics 
represented, are detailed in Section [3] alongside the fitting 
statistics and fitting process. Details of the calibration of 
the models and simulated observables are explained in full 
in Appendix [Bj where the online interface and tools are also 
described. The results of the fitting are presented in Section 
|4| where we also discuss individual objects and comparison 
with other works. Finally we conclude in Section [6] 



2 THE DATA 

We have selected and fitted a sample of LMS and BDs from 



Taurus presented in Guieu et al. (20071. Guieu et al. (20071 



present photometry for these objects and split them into 
those stars with infrared excesses (and by inference discs) 
and those without. For the purpose of this study we have 
used the same division. We have also performed a literature 
search for photometry and derived parameters of these ob- 
jects. This has allowed us to develop an understanding of 
how to model the uncertainties of our models during fitting 
and to explore temporal variations in the observations. 

Table [T] shows the short name used in this work and 
the 2MASS identification (from which the coordinates are 
easily found) to facilitate cross matching. Table p] shows the 
adopted photometry for each of our target objects and their 
sources. 

Table |2| also lists our estimates of the uncertainties 
for each observation. The uncertainties are not simply the 
values presented alongside the observations but also in- 
clude components which attempt to account for the dy- 
namical changes, (i.e. the fact the observations are non- 
contemporaneous) and inaccuracies in photometric trans- 
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Short Name 


2MASS ID 


Sources without infrared excess 


KPNO-Tau 4 


J04272799-I-2612052 


CFHT-Tau 15 


J04274538-I-2357243 


KPNO-Tau 5 


J04294568-I-2630468 


CFHT-Tau 16 


J04302365-I-2359129 


CFHT-Tau 13 


J04312669-I-2703188 


CFHT-Tau 7 


J04321786-I-2422149 


CFHT-Tau 5 


J04325026-I-2422115 


CFHT-Tau 11 


J04350850-(-2311398 


KPNO-Tau 9 


J04355143-f 2249119 


CFHT-Tau 2 


J04361038-I-2259560 


CFHT-Tau 3 


J04363893-I-2258119 


ITG 2 


J04380083-I-2558572 


Sources with infrared excess 


CFHT-Tau 9 


J04242646-I-2649503 


KPNO-Tau 6 


J04300724-I-2608207 


KPNO-Tau 7 


J04305718-I-2556394 


CFHT-Tau 12 


J04330945-I-2246487 


BDD399 


J04381486-(-2611399 


GMTau 


J04382134-I-2609137 


CFHT-Tau 6 


J04390396-I-2544264 


CFHT-Tau 4 


J04394748-(-2601407 


CFHT-Tau 8 


J0441 10784-25551 16 


BDD304 


J04414825-I-2534304 


BDD164 


J04442713-I-2512164 



Table 1. Table showing the short name used in this paper and the 
2MASS ID from which the coordinates can be found. The division 
of sources with and without infrared excess is from IGuieu et al.l 
l |2007| . 



formations. We note that although we have hsted the uncer- 
tainties in Table [2] we adopt a Bayesian approach whereby 
the observations are precise and the uncertainties are a prop- 
erty of the model. This is equivalent to asking which under- 
lying model, when combined with uncertainties representing 
our incomplete understanding of this model, will reproduce 
the precise observations. The construction of the uncertain- 
ties and how they are applied is described in detail in Ap- 
pendix |A] 



3 THE MODELS AND FITTING 

In order to model a combined star and disc system one re- 
quires a model of the stellar properties (luminosity, radius 
etc), the SED of the stellar photosphere and an initial model 
of the disc (structure, composition etc). The disc model is 
described as an initial model as some properties such as 
the vertical scale height or inner edge location may change 
(due to modelled physical processes) to some equilibrium 
state. Most existing studies of BDD systems involve param- 
eterisation in order to remove the requirement for detailed 
physical modelling. For instance commonly made assump- 
tions are that the central star emits as a black body, or as a 
naked photosphere without additional accretion flux. More 
subtly, one can simply pick a range of atmospheres at dif- 
ferent effective (Tch) temperatures, surface gravities (log((;)) 
and bolometric luminosities (Lboi) to create a fitting grid, as 
opposed to using stellar interior models to derive these prop- 



erties. Additionally, the inner edge and vertical structure of 
the disc can simply be prescribed analytically without in- 
clusion of the detailed physics of magnetic truncation, dust 
sublimation and vertical hydrostatic equilibrium. 



3.1 Dust Distribution 



We have adopted a similar dust model to Wood et al. ( 2002 1 



with the size distribution of dust particle given by 
n{a)da — Cia~'' X 63;^'"'"'°"' 'da 



(1) 



where n{a)da is the number of particles of size a (within the 
interval da), Ac is the characteristic particle size, with p and 
q simply used to control the shape of the distribution, d 
controls the relative abundance of each constituent species 
(j) in the dust. We have adopted the following values for 
each parameter q =3.5, p =0.6 with ttc =50 /im, also with 
an associated maximum and minimum grain size of amin =5 
nm and flmax =lmm. We set d =1 and adjusted the species 
using a grain fractional abundance (an equivalent process to 
that of Wood et al.|2002 l, using the solar abundances of As- 
plund et al. (2006 1. The subsequent opacities are constructed 



from tabulated opacities for each species. More detail of the 



dust model can be found in Mayne & Harries (20101 



3.2 The Grids 

In order to construct the stellar photosphere we interpolate 
the stellar interior models for a given age and mass and sub- 
sequently interpolate the stellar atmosphere models to the 
required Tcff, Lboi and \og{g). This method is in contrast to 
the traditional method of selecting spectra over a range of 
Teff, log((;) values and scaling the Lboi to match the observed 
fluxes a t wavelengths shortward of the peak in disc emission 

The accretion flux is then 



(~3Aim, |Dullemond et al.|2001 



added, as prescribed in |Mayne fc Harries (2010) , upon adop- 
tion of a rotation period (to derive the co-rotation radius). 
A disc component is then added and the radiative transfer 
through the system modelled to produce a simulated SED. 
The modelling of the BDD systems was performed using 
the TORUS radiative transfer code. TORUS is detailed in lHar-l 
ries| (|2000[) and was subsequently updated by [Harries et al.| 



(20041 and Kurosawa et al. (20041. torus uses the method 



of Lucy ( 1999 1 to solve radiative equilibrium on an AMR 



mesh. The code can also self-consistently solve the equation 
of vertical hydrostatic equilibrium and dust sublimation for 
the disc (iTannirkulam et al.l 120071). The TORUS code has 



been extensively benchmarked (for radiative transfer in discs 
against Pinte et al.|2009) . The addition of variable dust sub- 
limation is described in Mayne & Harries (20101. We have 
constructed two grids using different assumptions about the 
properties and physics within the star disc system. 



3.2.1 Ad-hoc Grid 

The first model grid includes the effects of accretion. We 
place the inner edge of the gaseous and dust disc at 
the co-rotation radius and prescribed an initial analytical 

^ In effect this is a 'by-eye' pre-fitting to select an input stellar 
photosphere. 
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vertical disc structure. The initial disc structure is con- 
structed using the analytical relationships (for the disc den- 
sity and scaleheight) from Mayne & Harries (20101, p — 
PofS-exp{-^[z/h{r)f) ex r~" and h = ho{r/Rl)'^ where r 
and z are the radial and vertical cylindrical polar coordi- 
nates. The surface density (E) profile from these equations 
is E(r) oc r^^~°'' , and /3 and a are set such that S(r) oc r^^ . 
This initial disc configuration follows the analytically pre- 
scribed disc structures of previous works, where SED fitting 
has been used to derive, for instance the /3 parameters or 
inner edge locations. These initial models are subsequently 
treated in one of two ways. Either, the vertical structure 
is fixed in the initial configuration and dust sublimation is 
modelled explicitly alongside a calculation of radiative equi- 
librium. Alternatively, hydrostatic equilibrium is enforced 
(via adjustment of the vertical structure and therefore /3) in 
concert with the radiative equlibrium and dust sublimation 
calculations. This grid is termed the ad-hoc grid and is an 
extension of the grid featured in Mayne & Harries ( 2010 1. It 



covers only BD mass objects for the ages of 1 and lOMyrs. 



3.2.2 Semi-empirical Grid 

The second grid is designed for non-accreting systems. In 
this CEise the inner edge is placed at the semi-empirically 
derived dust sublimation radius (only applicable for negligi- 
ble accretion rates, 



Whitney et al.^^2004J and the verti- 



cal structure analytically prescribed, as for the ad-hoc grid. 
In this grid, the vertical structure and inner edge location 
do not change as dust sublimation and vertical hydrostatic 
equilibrium are not included. This grid is termed the semi- 
empirical grid and has an expanded mass and age range 
covering masses from 0.01 to 1.40 M© and ages from 1 to 
10 Myrs. 

Note that our construction of the initial disc structures 
assumes that the dust is evenly distributed through the disc 
and follows the gas structure. We therefore implicitly ex- 
clude the possibility of dust settling in the disc for those 
models within the ad-hoc grid which enforce a hydrostatic 
structure on the gas distributions. However, the remaining 
models (i.e. a subset of the ad-hoc grid and all of the semi- 
empirical models) where the disc scaleheight structure is not 
required to fulfil hydrostatic equilibrium, can mimic the set- 
tling of the dusty component of the disc towards the mid- 
plane. Further complexity may be introduced by allowing 
the dust scaleheight to vary as a function of grain size (e.g. 
Tannirkulam et al.|2007 2008 1 although we believe that ad- 
ditional observable constraints (such as scattered light imag- 
ing see [Wood et ai:i[2002l [Watson fc Stapelfeldt||2007| for 
example) over and above SED fitting would be necessary 
to determine the degree of sedimentation. However, if dust 
settling or sedimentation is likely to be an important factor 
for a particular object then the semi-empirical grid or non- 
hydrostatic ad-hoc models are to be preferred. The details of 
which physical process are included in each grid are shown 
in Table IIE] 



^ This information is also available online at http://bd- server, 
astro. ex. ac .uk/calibration. 



3.2.3 Photometry 

In most cases parameters are derived from photometric ob- 
servations, either magnitudes or monochromatic fluxes, or 
conversions of the former to the latter. Therefore our SEDs 
need to be converted to apparent SEDs (via the applica- 
tion of a distance and extinction) and photometric observa- 
tions simulated. We have simulated such photometric obser- 
vations for a range of commonly used filter systems and the 
conversions, derivations of photometric measurements and 
definitions of the photometric systems are detailed in Ap- 
pendix IB] This information is important for observers who 
intend to use the fitting tool, as errors when adopting the 
incorrect photometric systems can be large (as mentioned in 
Section El and explained in Appendix \A\ . 

3.3 Fitting Statistic 

As shown in Table[2]we have collected a large range of photo- 
metric observations in both magnitude and monochromatic 
fiux systems. In order to perform a fit we must use consistent 
units and therefore we convert magnitudes to monochro- 
matic fluxes. This conversion is performed in a rigorous man- 
ner by reversing the process described in Appendix [B] used 
to derive a simulated magnitude and then applying the pro- 
cess required to construct monochromatic fluxes (again de- 
scribed in Appendix [B|. This allows us to be confldent that 
the models and data are being compared in a consistent 
manner. 

Practically, fitting is performed by maximising the like- 
lihood of the data being drawn from a given model. As the 
probability of the data being drawn from the ith model, Pm 



n^ 



Uqc 



rrij ^ 



o-J 



(2) 



where Xi and rrii are the ith data and model points (these are 
fiuxes from the SED, monochromatic fiuxes or magnitudes 
converted to fiuxes), and ai the uncertainty of the ith data 
point. The second part of the equation derives the likelihood 
from the observed upper limits xj is the jih upper limit 
with uncertainties aj , Mj the corresponding model fiux and 
where Q is given by. 



Q{x) 



-erf(^) 
2 ^^/2^ 



(3) 



where erf (x) is the error function. This in effect assumes that 
the a is taken from a measurement of zero, but allows neg- 
ative measured fiuxes. Ideally, one would provide the mea- 
sured flux and uncertainty as oppose to an upper limit, as 
this ensures the fltting routine does not have to make further 
assumptions about the data. In practice we choose to min- 
imise the s tatisti c — 21n(Pm), which as described in Naylor 
& Jeffries ( 2006 1 reduces to x^ for the case of data points 



with Gaussian uncertainties in one dimension (see also Cash 



1979 1. In practice therefore we minimise the following statis- 



tic, 



Apseudo 



-2\nil[e 



X + 



Uqc 



TTlj 



crj 



-21n(Q(- 



mj 



o-j 



(4) 



(5) 
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Grid Reference 



ad-hoc Grid 



semi-empirical Grid 



Physical Processes 



Stellar Interior 

Stellar Atmosphere 

Accretion Luminosity 

Vertical Structure 

Inner Dust Boundary 
Inner Gas Boundary 
Dust Scattering 



Dusty-'OO(l) 

AMES-Dusty(l) 

Yes 

VHE 

Static 

Variable Sublimation 

Co-rotation Radius 

Isotropic 



M, >=O.O2M0:BCAH98(2) 

M« < O.O2M0: Dusty-'00(i' 

Teff <1400K:AMES-Cond(3) 

1400 K< Teff <2500K:AMES-Dusty(i) 

Teff >2500K:NextGen(4) 

Yes-negligible 

Static 

Empirical Sublimation Radius 
Co-rotation Radius 
Full Mie Phase Matrix 



Parameters Covered 



Age (Gyrs) 

Mass (Mq) 

Teff 

log{5) 

Luminosity 

Accretion Rate log M 

Areal Coverage (%) 

Disc Mass (M,) 

Disc Outer Radius (AU) 

a in 

Scaleheight at lOOAU (AU)*^) 
Inclination (°) 



0.001 &; 0.01 

0.01-0.08 (A=0.01) 

Derived 

Derived 

Derived 

-12- -6 (A=-l) 

1, 10 



0.001-0.01 (A=0.001) 

0.01-0.10 (A=0.01), 0.10-1.40 (A=0.1) 

Derived 

Derived 

Derived 

-12 

10 



0.0, 0.01, 0.001 
100, 300 50, 100 

Derived (VHE), 2.00 (1.00), 2.10 (1.10), 2.25 (1.25) 2.00 (1.00), 2.10 (1.10), 2.25 (1.25) 
25 10 

0, 27, 39, 48, 56, 64, 71, 77, 84, 90 



References. (1) Chabrier et al. 1 2000 i, (2) Baraffe et al. (1998i, (3) AUard et al. 1 2000 i and (4) Hauschildt et al. 1 1999 i 



Table 3. Table describing the physical processes modelled in the two grids and the range of parameters the grids cover. (1) See explanation 
in Section^] (2) this initial scaleheight is allowed to adjust if vertical hydrostatic equilibrium is enforced. 



The online fitting tool then returns the A'^ best fitting mod- 
els to the user, with associated plots of Xpsoudo for each pa- 
rameter. Uncertainties are derived using the 68% confidence 
interval by converting the given x^ values to a probability 



(e- 



-(- 



'-)\ 



and integrating from the most likely value until 



68% of the total probability has been reached. Additionally, 
a Spearman's rank correlation coefficient is derived for the 
best fitting 10% of models for each pair of free parameters, 
giving an indication of degeneracy. If two parameters are 
found to have evidence of degeneracy 2D Xpscudo niaps are 
created and presented alongside the fits, although this is 
only an indication of degeneracy. 

The structure of the online fitting tools and com- 
ments on its use are detailed in Appendix IB] Appendix [B] 
also includes details of the other tools available at http: 
//bd- server .astro.ex.ac.uk/ such as an SED browsing 
tool and isochrone server. 



4 RESULTS 

In traditional SED fitting values of extinction and distance 
are set to single values (i.e without uncertainties) for a given 
fit, allowing the number of free parameters to be reduced 
by two. However, this discards useful information as the 
distance and extinction uncertainties are well founded con- 
straints. The depth (i.e. range of distances in a given pop- 
ulation) of a given group of stars contributes little to shifts 



within the HR-Diagram ( Mayne et al. 2007 Mayne & Naylor 



20081. However, when fitting monochromatic fluxes, or SED 
sections, to model SEDs variations in distance to individual 
objects can affect the flux significantly. For instance when 
fitting objects within Taurus the distance is usually set to 
140 pc, or in some cases a modest distance uncertainty of 
zblOpc is used (Bertout et al. 111999 1. However, the distance 



range, from the most precise derivations of distances to three 
member stars, is 161±0.9, 147 and 130 pc, derived from in- 
terferometric observations ( Torres et al.|[2009 l. Therefore a 
more realistic distance range is 130-165 pc, which can lead to 
fiux variations of up to a factor of 1.6. Assuming a distance of 
140 pc for all objects in the Taurus region allows us to reduce 
the number of free parameters, but is less desirable than us- 
ing a realistic range of distances. Uncertainties should also 
be adopted to apply to the extinction to each target ob- 
ject. The extinctions shown in Table[2]have all been sourced 
from Guieu et al. (20071 (for consistency), where they are 



adopted as single values (without uncertainty). These val- 
ues, in turn, are derived in Guieu et al. (20061 where an 



uncertainty of ±0.8 mags is quoted, which we have adopted. 
If we do not include such strongly physically motivated ex- 
ternal constraints, like distance ranges, it is likely that, due 
to degeneracies of the parameters, the remaining free pa- 
rameters will be optimistically constrained. 

Figure [l] shows the infrared and near-infrared colour- 
magnitude (left panel) and colour-colour (right panel) dia- 
grams (CMD & CoCoD) for targets from Table |2] A 1 Myr 
isochrone constructed using the interior models and atmo- 
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Harries (20101 showed that the effects of accretion and disc 
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Figure 1. Figure showing the stars with (asterisks) and without 
(empty circles) infrared excesses from Table ^ The left panel 
shows a J,H - K CMD and the right panel a ([3.6]-[4.5]), 
([5.8] — [8.0]) CoCoD. The dashed line on the left panel shows an 
isochrone at 1 Myr for masses between 0.01 and O.OSMq adjusted 
to a distance of 140 pc and an extinction of Ay = 1- 



spheric models of |Chabrier et ah] ( |2000[ ) adjusted to a dis- 
tance of 140 pc and an Ay = 1, is overlaid on the CMD (left 
panel) . 

The left panel of Figure [T] shows that the stars with (as- 
terisks) and without (empty circles) discs lie redward and 
brighter than the illustrative isochrone. For the stars with- 
out discs several factors lead to a shift away from the model 
isochrone. Firstly, some objects lie behind extinction greater 
than the illustrative value selected (i.e. Ay >1.0). Addition- 
ally, some stars without discs have luminosities above that 
predicted by the 1 Myr isochrone at a given colour or H — K 
(for instance ITG 2, discussed in more detail in Section 
4.2.1 1. Finally, some stars will also have a mass greater than 
the maximum mass included in the illustrative isochrone 
(O.O8M0, for instance CFHT-Tau 5 also discussed in more 
detail in Section 4.2.1 1. Of course, for the stars with discs the 



same shifts will be present in addition to reddening of the 
system caused by the reprocessing of stellar light by the disc. 
As expected, the stars with infrared excesses are also clearly 
distinguishable in the right hand IRAC CoCoD, and the 
delineation agrees with the predicted values of ([3.6] — [4.5]) 
> -h0.21 and ([5.8]-[8.0]) > -h0.13 ( |Mayne fc Harries|2010| ). 
It is clear that the stars in Figure 111 do not lie along the 
isochrone, and therefore are poorly represented, in a CMD, 
by a simple coeval (naked photosphere) population at 1 Myr 
at a single distance and extinction with negligible accretion. 
The meaning and veracity of the inferred "age spreads" has 
been discussed at length in the literature (for a summary see 
Mayne fc Naylor|2008 1 . Several studies of classical T-Tauri- 
stars (CTTS) have shown that active accretion and varia- 
tions in distance cannot fully explain the spreads of stellar 



populations (Burningham et al. 2005 Da Rio et al. 20101. 
However, accretion, both past and current, has been iden- 
tified as a major contributor to spreads in both luminosity 
and radii of young stars, and its effects on the magnitudes 
and colours of brown dwarfs with discs has been shown to 



presence on characterisation are more significant for lower 
mass systems. 

A full discussion of the veracity of age spreads is beyond 
the scope of this paper. However, it is clear that fitting for 
an assumed coeval population can be performed by adopting 
a range of "isochronal ages" , whether they represent a stars 
true age or not. Additionally, a range of accretion rates and 
inclinations (leading to differing levels of occulation of the 
star by the disc) may be sufficient to provide the spreads 
observed, with the addition of a small intrinsic age spread 



which is less than the disc dissipation time (Littlefair et al. 
20TH [Jeffries et al.|201l"| . 



4.1 Literature Stellar Parameters 

For each of the objects within our sample we have attempted 
to find all the stellar parameters that authors have either de- 
rived or adopted for modelling and analysis of these objects. 
Table HI presents, for each object, each value of the stellar 
parameter we have found, with the reference in superscript 
parenthesis (a legend is included in the table caption). 

Table|4]does not contain information on how each quan- 
tity was derived (or indeed whether it was simply adopted), 
which is of course important. However, in cases where pa- 
rameters were simply adopted in order to calculate or derive 
a dependent parameter it is still important to note the val- 
ues of the input, which will, of course, determine the value 



of the derived property. For example Muzerolle et al. ( 2005 1 



adopted approximate values for stellar parameters such as 
Tcft and mass, in order to model accretion and derive accre- 
tion rates for such objects, whereas, Teff for these objects 
were derived from fitting to spectral type templates in |Luh-| 
man ( 2004 1. Of particular note is that almost all of the stel- 
lar masses derived for these objects rely, in some way, on 
isochrones, which are known to be unreliable for low mass 
stars (see discussion in Mayne & Naylor 20081, under the 



assumption of coevality at an age of 1 Myr, and crucially, 
often at a set distance of 140 pc. 



4.2 Naked Stars 

The naked stars, which we assume are not accreting (sup- 
porting evidence is presented in Table [7| discussed in Section 
4.4.1 1, were fit using the semi-empirical grid to allow a range 

There 



be significant (Mayne & Harries 20101. Indeed Mayne & 



of isochronal ages as observed in the CMD, Figure! 

is no difference, except the accretion luminosity, between the 

naked photospheres used in the ad-hoc and semi-empirical 

grids, although the latter includes a wider range of ages and 

masseq_] 

Figure 12] shows an example of a fit to CFHT-Tau 16 us- 
ing the semi-empirical grid of models and the observations 
presented in Table [2] We have fitted this star in three differ- 
ent ways. Firstly, we have fitted the star allowing all of the 



^ We note that these values should not be adopted as the real 

ages of the system. As discussed in Section |4] ages derived from 

isochrones are currently inaccurate for LMS and BDs. 

■* Of course, the models of systems at higher masses will include 

stellar atmospheres and interiors appropriate for the higher T^g 

values. 
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Figure 2. SEDs of the best fitting models from semi-empirical 
grid, showing fits to the entire grid (dashed green line), using a 
set distance and extinction (dashed black line) and a restricted 
grid (solid black line). The monochromatic fluxes derived from 
the best fitting model from the restricted grid are shown as red 
circles and the observations as blue crosses (fluxes) or triangles 
(upper limits) with uncertainties shown as bars. Note, the fits for 
this object are closely matched and therefore the resulting SEDs 
show little deviation, meaning the lines for each fit are essentially 
over plotted and some may not be visible. 



parameters of the model grid to be free. Secondly, we have 
fitted the star by setting the distance to 140 pc and extinc- 
tion to the single values from Tablep|(following |Guieu et al.| 
2007). Finally, we have performed a fit using observations, 



shown in Table HI to constrain the model grid parameters 
alongside distance and extinction uncertainties, fn the case 
of CFHT-Tau 16, all of the fits yielded similar results. The 
resulting x^ parameter space for the case where the param- 
eters are all free is shown in Figure |3] 

Figure [2] shows that in the case of CFHT-Tau 16 the 
best fitting SEDs, for each fit type, are well matched. This 
is reflected in the consistency of the derived parameters, 
for this object, shown in Table p] (described later). How- 
ever, Figure [3] shows that even in this excellent case some 
parameters remain formally unconstrained when the confi- 
dence intervals are examined. The 68% confidence intervals 
for the extinction and mass parameters are within the fit- 
ted parameter space, but the the distance upper bound is 
unconstrained. The distance ranges selected, however, are 
strongly constrained by the range of distances found to ex- 
isting members and should not be extended. It is usually 
inadvisable for an object to be fitted to the entire grid, since 
degeneracies in the parameter space can lead to derivations 
of incorrect parameters. For instance the effects of increas- 
ing mass, decreasing distance and increasing accretion rate 
all act to increase the luminosity of the system and therefore 
can lead to compensating errors, especially when combined 
with adjustments in extinctions. This is especially signifi- 
cant for lower temperature stars as the region of the SED 



sensitive to stellar and disc properties overlap ( Clarke fc 
Pringle||2006[ ). Indeed payne fc Harries] ( |2010[ ) show that a 
population of accreting BDD systems can appear both spec- 
troscopically and photometrically similar to a population of 
CTTS under heavy extinction and at a different distance. 



For the remaining sample of naked BDs we have fitted 
using the final two methods applied to CFHT-Tau 16, i.e. 
fitting using a set distance and extinction, and fitting using 
a restricted grid with all parameters constrained by obser- 
vations (in other words we do not fit using a grid where 
all the parameters are free). Figure |4] shows the best fitting 
SEDs for each of the naked objects. The solid and dashed 
lines show the fits allowing distance and extinction to vary 
over the prescribed ranges and setting them to fixed val- 
ues, respectively. The model fluxes, for the restricted case, 
are shown as circles and the observations as crosses, with 
triangles denoting upper limits. 

The parameters derived for the restricted fltting grid 
are shown in Table [5] Table [5] shows the best fitting values 
for the independent variables with 68% confidence intervals 
and the dependent variables associated with the best fitting 
model. Some of the confidence intervals indicate a 0.0 in- 
crement in either direction, this is sometimes caused by the 
best fitting value reaching a boundary (for example 165 pc in 
the range 130-165 pc as is the case for CFHT-Tau 16). This 
may indicate that the parameter is poorly constrained, how- 
ever, as the boundaries for the fitting grid have been chosen 
using external literature constraints this is not necessarily a 
problem. 

Comparison between the values derived in Table Island 
the literature values (Table p| show reasonable agreement. 
For most objects our derived masses match, within the un- 
certainties, those found in the literature and the effective 
temperatures match to within ~100K, with the exception 
of ITG 2, CFHT-Tau 5 (which are discussed in Section [4!2T| 
and KPNO-Tau 4. For KPNO-Tau 4 the discrepancy in the 
effective temperatures is large with a difference of ~350 K. 
This is most likely caused by this object existing at the 



boundary in our grid between the AMES-Dusty (Chabrier 



et al.||2000[ ) and A MES-Cond ([Allard et al.||2000[ ) models, 
which as shown in |Mayne fc Harries] ( |2010[ ) leads to a sig- 
nificant disconuity in photometric magnitudes. 

For most cases allowing the distance to change over the 
observed range, and adding an uncertainty to the extinc- 
tion, results in parameters which are consistent (within the 
uncertainties) with those derived with a set distance and ex- 
tinction (excepting ITG 2 and CFHT-Tau 5). However, the 
best fitting values are often slightly different, meaning the 
derived dependent variables do vary. For example changes 
in the T^s are, for the most extreme cases of CFHT-Tau 11 
and KPNO-Tau 9, 134 or 289 K respectively, however for the 
remaining fits they are below 100 K. 

4.2.1 Problem Fits 

We were unable to achieve acceptable fits using a restricted 
parameter space, to the objects ITG 2 and CFHT-Tau 5. 
For ITG 2 setting the distance and extinction and allowing 
the rest of the grid parameters to run free also resulted in a 
poor best fit model. For CFHT-Tau 5, however, an accept- 
able fit was found using a set distance and extinction but 
allowing the stellar mass to increase significantly. To explore 
the problems with these objects we have also fitted them us- 
ing the entire grid. The resulting best fits are presented in 
Figures [5] and |6] alongside those already presented in Figure 
[4] For ITG 2 the only way we could derive a best fitting SED 
close to the observed data points was to allow the mass of 
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Figure 3. Figure showing the x^ distributions for the fit to the whole semi-empirical grid for CFHT-Tau 16, where the dashed lines 
show the 68% confidence limits. For each plot the remaining parameters have been held at their best fitting values. 
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Figure 5. SEDs and monochromatic fluxes for ITG 2 where sym- 
bols have the same meanings as Figure [2] 
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Figure 6. SEDs and monochromatic fluxes for CFHT-Tau 5 
where symbols have the same meanings as Figure [2] 



the star (M«) to increase to 0.2 M© and allow the extinction 
and distance to increase to ~5 Av and 165 pc respectively. 
The resulting star has a much higher effective temperature of 
3193 K, log(3) of 3.40, radius of 1.45 Rq and a luminosity of 
log (L/Lq)=— 0.700. These values are inconsistent with the 
literature properties presented in Table |4] For CFHT-Tau 5 
the best fitting SED returned values of, M«=0.1 M©, dis- 
tance=140.50pc, Av=0.50, log (L/L0)=-1.16O, Teff=3002, 
log(5)=3.45 and a radius of O.968R0. Again, as with ITG 2, 
these parameters are inconsistent with the literature values 
presented in Table [4] 

These objects are at much higher luminosities, for their 



given masses, than the youngest (i.e. 1 Myr) isochrone. This 
means that these stars are poorly fit by both of our grids. 
Hence, we can only achieve these luminosities by increasing 
the stellar mass and compensating for the SED shape by 
increasing the extinction. Table [4] shows that CFHT-Tau 5 
has a similar luminosity to CFHT-Tau 7, however, the latter 
object has a significantly higher mass and is, therefore, not 
affected by this problem. This limitation means that for very, 
"isochronally" young stars we will be unable to reproduce 
their SEDs with our current grid. We are currently running 
a grid of models which are divorced from isochronal theory 
and will extend to this region of parameter space (explained 
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Figure 4. The best fitting SEDs for tlie naked objects, all from the semi-empirical grid. The solid black lines are the SEDs from the 
restricted fits and dashed lines show the best fitting SEDs when the distance and extinction have been set. The monochromatic fluxes 
derived for the best fitting restricted model are shown as circles and the observed fluxes, (from Table[2]l are shown as crosses, or triangles 
for upper limits, with the bars denoting the uncertainties. It can be seen that we failed to find a satisfactory fit for ITG 2 and CFHT-Tau 
5 denoted by (*). 



in more detail in Section [m. It is interesting to note that 
the SEDs of such objects cannot be reproduced if we allow 
the accretion rates to increase, even though they are naked 
stars (some supposed naked stars have been found with non- 
negligible accretion rates|Kennedy fc Kenyon|2009[ these are 
discussed in Section 4.4.1 1. The accretion diagnostics for the 
sample have been collected (Table[7| and are discussed later 
in Section 144. II 



4.3 Stars with Discs 

For the BDD objects we have fitted to the ad-hoc grid to 
allow for a range of accretion rates, again using either a set 
distance and extinction or the restricted grid with varying 
distance and extinction. One important note is that, as dis- 
cussed in Section r3.2.1| the observed dust disc stuctures may 
not appear in hydrostatic equilibrium due to grain settling. 
Therefore, we have included the ad-hoc models with pre- 



scribed disc structures in the fitting process. This may also 
be a reason to favour the semi-empirical models for systems 
with significant grain settling, but with negligible accretion 
rates. 

Figure [7] shows the best fitting SEDs for the BDD ob- 
jects from the ad-hoc grid. As for Figure |4] the solid and 
dashed lines show the fits allowing distance and extinction 
to vary over the prescribed ranges (the restricted fits) and 
those with distance and extinction set at fixed values, respec- 
tively. The model fluxes, for the restricted case, are shown 
as circles and the observations as crosses, with triangles de- 
noting upper limits. 

Figure [7] shows that the fits do not change significantly 
upon fixing the distance and extinction, with the exception 
of CFHT-Tau 8. It is clear that all observations are well fit 
with accreting 1 Myr disc structures drawn from the ad-hoc 
grid, again excepting CFHT-Tau 12. The best fitting param- 
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Figure 7. The best fitting SEDs for the BDD objects, from the ad-hoc grid. The soUd black fines are the SEDs from the restricted fits 
and dashed fines show the best fitting SEDs when the distance and extinction are fixed. The monochromatic fluxes derived for the best 
fitting restricted model are shown as circles and the observed fluxes, from Table [2] are shown as crosses, or triangles for upper limits, 
with the bars denoting the uncertainties. It can be seen that we failed to find a satisfactory fit CFHT-Tau 12 and, to a lesser extent, 
CFHT-Tau 4, which are denoted by (*). 



eters derived for the BDD systems using a restricted grid 
drawn from the ad-hoc model grid are presented in Table [6] 

We found that fixing the distance and extinction only 
affects the derived best fit parameters significantly in a few 
examples when using the ad-hoc grid. For example, the best 
fitting mass for BDD304 increases by O.O2M0 when the dis- 
tance and extinction are set. For most objects the derived 
accretion rates are consistent with those found in previ- 
ous studies (see Table ml . Additionally, in all cases except 
CFHT-Tau 12 and CFHT-Tau 9 the parameters derived for 
the targets using the accretion models (i.e. ad-hoc) are con- 
sistent with the spread of values presented in Table HI It is 
clear from Table [6] that we are unable to differentiate be- 
tween accretion rates ranging from logM=— 12 to —9. This 
is due to unavailability, or poor quality of, optical photome- 
try which covers the shorter wavelength range where accre- 
tion flux is the dominant contributor to the SED. The de- 
rived accretion rates are discussed in detail in Section |4AT] 



We also fitted the BDD objects to the semi-empirical 
grid, where accretion rates are negligible but a larger range 
of ages and masses is available. For all BDDs, excepting 
BDD399, BDD304, BDD164 and CFHT-Tau 4, acceptable 
fits could be found using the non-accreting models. This is 
as the effects of a disc and accretion on observations of a 
BDD system can mimic a star of higher mass or younger 
age ( Mayne fc Harries|2010 \ . Therefore observations of sys- 
tems with discs, can be fit using SEDs of accreting and non- 
accreting models. For instance for CFHT-Tau 8 and KPNO- 
Tau 6 best fitting masses of 0.1 and 0.05 M© (and ages of 
5 Myrs) were found with the semi-empirical grid, compared 
to 0.08 and 0.03 found using the ad-hoc grid, and 0.078 and 
~0.025 M0 (for an age of 1 Myr) found in the literature (see 
Table HI . It is perhaps surprising that even in these cases, 
where a large range of photometric measurements are used, 
that there can still be some confusion as to whether a source 
is accreting, or if it is in fact a more massive and younger ob- 
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ject. This suggests, as detailed in Mayne & Harries (20101 



that selection of stars from CMDs and subsequent inves- 
tigation into their mass and accretion rate relationship is 
unreliable. In the case of CFHT-Tau 12 the fit is actually 
better when using the semi-empirical grid, however several 
models cannot be fit without allowing for accretion flux, i.e. 
BDD399, BDD304, BDD164, GM Tau and to a lesser extent 
CFHT-Tau 4. 



4.4 Derived Parameters 

The main differences between our models of BDD systems 
and those used in previous studies is the additional mod- 
elling of accretion luminosity, vertical hydrostatic equilib- 
rium and dust sublimation. The main derived parameters 
which indicate the veracity of these mechanisms, when com- 
pared to observations, are the accretion rate, inner edge lo- 
cation and flaring parameter. Without the inclusion of these 
physical mechanisms it is clear from previous works and from 
the fits in this paper that reasonable fits can still be achieved 
and parameters derived. However, if important physical pro- 
cesses are neglected it is questionable whether these parame- 
ters are reliable and whether we are learning anything about 
the physical mechanisms governing these systems. 



4.4-1 Accretion 

Accretion rates of logM< —9 contribute little to shifts 
within the HR-Diagram ( Mayne fc Harries||2010 l. However, 
when fitting monochromatic fluxes or SED sections to model 
SEDs variations in accretion rate below this threshold can 
affect the flux significantly (Mayne fc HarrieS[ 2010). This 



can apply to objects with and without discs e.g. |Kennedy| 
& Kenyon (20091 find 43 stars in their sample of 1253 are 



actively accreting whilst no disc is detected. 

The adopted magnitudes (from Table |2| do not extend 
across the optical regime, i.e. magnitudes at U, B and V. 
This does mean that a large section of the SED remains un- 
constrained. Fitting has been performed by adopting mag- 
nitudes or deriving fluxes from photographic plate observa- 
tions, such as the NOMAD catalogue (Zacharia s et al.|2005 1 
by |Bouy et al.|p008[ ). However, the spectral responses of the 
photographic plates are poorly constrained and therefore 
subsequent transformation to another photometric system 
or to a monochromatic flux are unreliable (see Appendix 
B2L However, as accretion flux can contribute signiflcantly 
in the optical regime allowing this parameter to run free 
may lead to best fitting models having elevated or unrealistic 
levels of accretion. Therefore, we have collected photometry 



from Zacharias et al. (20051 for as many of our objects as 



possible alongside Ha equivalent widths and derived accre- 
tion rates, which are presented in Table [T] As many of our 
objects cannot be satisfactorily flt without accretion flux 
these values can be used to support their status as active 
accretors. These values can only be used as a guide to add 
support for an accreting or non-accreting model. The op- 
tical photometry is quite unreliable and it is very difficult 
to compare measured widths of Hq made with different in- 
struments or over different epochs due to their variability 
( Jayawardhana et al.||2006| , so one cannot use these values 
as strong constraint. 



Table [7] shows accretion diagnostics in the form of Hq 
equivalent widths for our targets as found by various au- 
thors, and derived accretion rates for some objects. Table [7] 
also shows our derived accretion rates for the BDD objects. 

We have assumed negligible accretion rates for all of 
the BD systems and Table [T] shows that this is a good as- 
sumption for these objects. For the BDD objects, in most 
cases, our accretion rates are consistent with those derived 
in the literature. Alternatively, where we have derived non- 
negligible accretion rates the EW(Ha) and blue photometry 
support this. Several individual cases, however, will be in- 
vestigated further in Section [4. 4. 3| 

4-i-2 Disc Structure 

The ad-hoc grid includes a state-of-the-art treatment of dust 
sublimation and vertical hydrostatic equilibrium. Therefore, 
one would expect to find significant differences between the 
inner edge locations, inclinations and fiaring parameters de- 
rived in this work, and those derived by studies using an 
analytical or semi-empirical approach (for instance placing 
the inner edge at a predefined sublimation radius). 

Table [8] shows our best fitting inclinations, inner edge 
locations (converted to AU) and /3 (1 — a, the flaring param- 
eter) compared against those derived in |Guieu et al. (20071 
for the BDD targets, where availablaj The disc structures 
and surface density profiles are varied in|Guieu et al. (20071 



but the corresponding value of the surface density exponent 
of the best fitting model is not published, so comparisons in 
TablelSlare incomplete. Our best fitting total disc masses are 
presented in Table [6| for each of the BDD systems. Whereas 
the models used in |Guieu et al. ( 2007 1 adopted a single disc 
mass of 1 Mj. 

Tablelslshows that, in general, our inner edge locations, 
or more precisely the position of our dust sublimation radii, 
are often much closer to the star than those found in lGuieril 



et al. ( 2007 1 . This is caused by the inclusion of accretion flux 



in our models and realistic dust sublimation as opposed to 
the prescription of a vertical wall. As discussed in |Mayne| 
& Harries (20101 the resulting inner edge structures are 



curved, consistent with detailed studies ( Tannirkulam et al 



2007 1 of the inner edges of circumstellar discs, which leads 



to significant changes in the way the SED around l-lO/xm 
changes with inclination. This means we are able to achieve 



fits aesthetically similar to those of Guieu et al. ( 2007 1 but 



with significantly different inner edge locations. It is how- 
ever difficult to perform a complete comparison as |Guieu| 
et al. (20071 do not publish their best fitting inclinations. 



Additionally, Guieu et al. ( 2007 1 only model inner edge po- 



sitions of 0.015-0.3 AU from the central star, which does not 
extend low enough to match some of our models. 

4--i-3 Specific Cases 

Further comparison and investigation into the structures of 
the discs can be performed in a few cases where either, 
the fitting process has revealed interesting results or further 
studies have been performed on the object in question. 



^ [Guieu et al.| | |2007| do not publish their best fitting inclinations 
or the fits for CFHT-Tau 12 and BDD399 
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Name 



Equivalent Width of H^ (A) 



log A/ 



B, 



plate V'platc log Af * 



Sources without infrared excess 



CFHT-Tau 16 -16.76(^) 

CFHT-Tau 15 -18. go'^' 

CFHT-Tau 13 -S.GO'i' 

CFHT-Tau 11 -45.07(1) 

CFHT-Tau 7* -8.63(i),5.5(*) 

CFHT-Tau 5 -29.84(i) 

CFHT-Tau 3 -10.5W,50±4 or 65±l(3), -55±4('^) <-12('*) 

CFHT-Tau 2 ll±l(8),-13±4C^) 

KPNO-Tau 9 20(2) 

KPNO-Tau 5 30(3) 

KPNO-Tau 4 150(2) -38.4W-158.08(i) <-12('*) 

ITG 2 



20.31 
18.83 



20.70 
20.56 



17.97 
17.96 



<-12 
<-12 
<-12 
<-12 
<-12 
<-12 
<-12 
<-12 
<-12 
<-12 
<-12 



Sources with infrared excess 



CFHT-Tau 12 
BDD399 
BDD304 
BDD164 
CFHT-Tau 9 
CFHT-Tau 8 
CFHT-Tau 6 
CFHT-Tau 4 
CM Tau 
KPNO-Tau 6 
KPNO-Tau 7 



-79.4(1) 

-47(4) 

-233.7(*) 

-9.95(1) 
-52(1) 

-102(''), -63.7(1) 
-129.3(''), 69±4(3) 

350(2.3) ^_207.9l(i), -41. l(3) 
300(2.3) -31. !('*) 



-10.8(-i) 

-11.3('i) 20.33 

-ll.OSC'), -9.7('') 19.27 

20.59 



-11.3('i) 
-11.3(*) 

-8.7(6) 

-11.4(4.5) 
-11.4(4) 



19.99 
18.39 



17.89 



17.45 



ir|-|-3.00 

1^-0.00 

0-1-0. 07 

''-1.00 

Q-l-1.03 

^-1.00 

0+0.55 

"-3.34 

Q+0.13 

^-3.00 

H.03 

-1.00 

a-l-0.06 

-3.00 

a-l-0.14 

-3.00 

3+0.63 

'-3. 11 

",+0.05 

-2.00 
-1+0.90 

-2.00 



-10"' 
-9^ 

-9^ 
-8^ 
-10^ 
-lOJ 



Table 7. Table of accretion indicators for sources with and without infrared excesses as defined in|Guieu et al.|(|2007j. (*) the log A'/in 
column 6 is from this work. References: ( l)|Guieu e t al. (2006), (2) Bri cefio et al.|||2002|, (3 ) Jayawardhana et al. (2003), (4) MuzcroUe] 
et al.|(|2005], (5)|Herczeg fc Hillenbrand|p008| l, (6) |White k. Basri,(,2003|l , (7) [Martin et al.| ( J2001J , (8) |Kraus fc Hillenbrand, ( ,2009 J and 
Bouy et al. (|2008|. Bpi^te and Vpiate photometry is fromlZacharias et al.l (|2005|. 



CFHT-Tau 12 is better fit using a semi-empirical model 
with a negligible accretion rate, rather than an accreting 
model from the ad-hoc grid. Figure IS] shows the fits for 
CFHT-Tau 12 using both grids. Figure [8] shows that the 
semi-empirical grid yields a much better fitting SED in all 
cases. Most parameters, and their uncertainties, are consis- 
tent between the fits in the two grids. However, the best 
fitting mass is marginally increased, the distance and incli- 
nation decreased, and the age moves from 1 Myr to 2 Myrs. 
Only by allowing the accretion rate to move to logM=— 7 
can we get an acceptable fit using the accreting, ad-hoc mod- 
els (dashed green line left panel, Figure [8|. Such an extreme 
accretion rate is unlikely but there is little additional in- 
formation in Table IT] to constrain this, other than a non- 
contemporaneous H(j equivalent width measurement. It is 
interesting to note that CFHT-Tau 12 has no measurements 
at wavelengths longer than 24 /xm (see Table [2| leading to 
poor constraint of the cooler dust disc at larger radii, how- 
ever our best fitting disc mass is Mdisc=0.001M«, with no 
acceptable fit to an O.OIM* disc model found. 



BDD399, BDD164, BDD304, CM Tau and CFHT-Tau 
4 are all poorly fit when using the semi-empirical grid, (i.e. 
assuming negligible accretion and allowing age to vary) . This 
is not unsurprising as for each of these objects there is ev- 
idence of non-negligible accretion rates, as shown in Table 
[7] If such objects are fitted without accretion luminosity 



the mass of the central star is increased and its distance 
decreased, to increase the short- wavelength luminosity. 

Figures |9] [lO) [TT] [12] and [13] show the fits from both 
grids (semi-empirical and ad-hoc) for all fit types (i.e. using 
the entire grids, refining the fit and finally setting only the 
distance and extinction). 

For BDD399 fitting to the entire semi-empirical grid, 
still does not result in a satisfactory fit. As stated in |Luh-| 
man et ah] ( [2007 1 the photosphere is too faint for any model 
of a young BD. However, as we can see in the left panel of 
Figure [9] the fit is much more consistent for all fit types in 
the ad-hoc grid, and indeed yields similar parameters. The 
derived accretion rate for this object is log-M=— SI^qq (for 
the restricted fit) which is greater than the literature result 
of —10.8 (MuzeroUe et al. 2005). This object was studied in 



detail in Luhman et al. (2007) where an inclination of 70° 
was found with an inner disc radius of 58R* and outer radius 
of R=20-40 AU. The central star of jLuhman et al.l ( [2007 l 

had a luminosity of —1.2218 log ( j^ ), a mass of O.OSMq 
and a radius of 1.O2R0. We find an inclination angle of 64°, 
an inner disc radius of 6.4R* and an outer disc radius of 
100 AU, with a corresponding star of luminosity —1.3400 
log (jf-\ mass of O.O8M0 and radius 0.829Ro. We derived 
a best fitting disc mass for this object of Mdisc=0.001M«. 
Whilst our inclinations are consistent with that of Luhmanl 



et al. ( 2007 1 our fit is hampered by the fact that the luminos- 
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Figure 8. SEDs and monochromatic fluxes for CFHT-Tau 12 where symbols have the same meanings as Figure [2] 
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Name 


Source 


Inclination (°) 


/3 (0=VHE)(i) 


Rsub (AU)(2) 


Sources with infrared excess 


CFHT-Tau 9 


This paper 


&itir 


(1.25) 
1.125 


0.0081 


Guieu et al. 


poor 


0.032-0.067 


KPNO-Tau 6 


This paper 


m+_ir 


(1.25) 
1.0-1.125 


0.0036 


Guieu et al. 


poor 


0.032 


KPNO-Tau 7 


This paper 


7ll?f 


1.25 (0) 
1.0-1.125 


0.0043 


Guieu et al. 


l|2007 


0.015 


CFHT-Tau 12(*) 


This paper 


6Atlf 


1.25 (0) 
Not Published 


0.0058 


Guieu et al. 


l|2007 




BDD399 


This paper 


g4+0.00 


1.10 
Not Published 


0.025 


Guieu et al. 


poor 




GM Tau 


This paper 


4g+8.01 


1.10 
1.0 


0.025 


Guieu et al. 


poor 


0.015-0.067 


CFHT-Tau 6 


This paper 


48l«f 


1.25 (0) 
1.0-1.125 


0.064 


Guieu et al. 


l|2007 


0.067-0.14 


CFHT-Tau 4 


This paper 


39t?2*^ 


1.10 
1.0 


0.025 


Guieu et al. 


l|2007 


0.14-0.3 


CFHT-Tau 8 


This paper 


7-1+0.76 
'-'--49.43 


(1.25) 
1.0-1.125 


0.025 


Guieu et al. 


poor 


0.015-0.067 


BDD304 


This paper 


7-1+2.17 
'-^-71 


1.25 (0) 
1.125 


0.022 


Guieu et al. 


poor 


0.015-0.032 


BDD164 


This paper 


„+56 
^-0 


1.10 
1.125 


0.025 


(|Guieu et al.|2007 


0.032-0.067 



Table 8. Table of disc structure parameters for sources with infrared excesses as defined in [Guieu et al. (20071. Here our Rsub radii have 
been converted to AU and a /3 parameter constructed using 1 — a (see Sectional. Inclinations from our work arc included but inclinations 
are not published in [Guieu et al.| l |2007| |. (1) As the /3 variables can take only four values uncertainties are meaningless we simply quote 
the best fitting values and if models with the alternate value were found to have x^ values within the 68% confidence interval their values 
are quoted in brackets. (2) the Rsub from our work is positioned according to the balance of stellar flux and disc temperature, the range 
chosen in [Guieu et al.[ l |2007[ l is 0.015-0.3 AU. 



ity of the stellar photosphere of BDD399 is under-luminous 
when compared to isochronal models (a grid is in prepara- 
tion to improve on this limitation). As Luhman et al. (20071 
obtained images of the disc directly their derived disc pa- 
rameters should be viewed as more reliable. 

For BDD164 a satisfactory fit can only be achieved with 
an accreting model, however in this case our derived accre- 
tion rate of logM=— 81334 is consistent with the values of 
-11.03 and -9.7 ( |Muzerolle et al.|[2005irBouy et al.||2008| ), 
albeit with the best fitting rate elevated over the literature 



![oo, compared 



values again. Bouy et al. (20081 studied this object in detail 



using multi-wavelength photometry (although as noted in 
Section [4.4. 1[ they used NOMAD photometry in the optical 
regime which is derived using photographic plates in this 
case). They found the best fit for a disc with inner radius 
at 0.04 AU and a /3 value of 1.1 or 1.2 (consistent with an 
a of 2.10 or 2.20 in this work) for a central star with mass 
O.O45M0. We find an a of 2.10 and an inner disc radius of 
0.025 AU for a star with mass 0.07-0.08Mq, with a best 
fitting total disc mass of Mdisc=0.01M». 

For BDD304, Figure [TT] again reveals a requirement of 
an accretion fiux when this object is fitted and comparison 
of the derived parameters from Tables [6] and [T] shows that we 



derive an elevated accretion rate, log Af=— 9 

to -11.3 ( |Muzerolle et al.|2005| ). We find a best fitting disc 



mass of Mdiac=0.001M* for this object. 

For GM Tau, even using the accreting models of the ad- 
hoc grid we are unable to achieve a completely satisfactory 
fit, as shown in Figure [121 This is caused by the fact that the 
mass of GM Tau is estimated to be at or above the higher 
mass boundary of the ad-hoc grid (M« — O.O8M0, compared 
to 0.078, see Table |4|. Indeed, this can be seen when fitting 
using the semi-empirical grid where the mass derived is much 
higher. It is interesting to note that if the monochromatic 
fiuxes were fitted to the entire grid (only setting the distance 
and extinction), as shown in Figure 12 a satisfactory fit can 



be obtained whereby the object looks very similar to a higher 
mass object at higher inclinations, Mt=O.8M0 at 64° (as 
highlighted by [Mayne fc Harries||2010| ). This highlights the 
problems of using fitting tools such as the one featured in 
this paper and that of [Robitaille et ah] ( |2006[ ) without pre- 
selecting, or refining the parameter space (this point is also 
emphasised in Robitaille et al.[[2006 |. The best fitting disc 
mass derived for this object was Mdisc=0.01M«. 

As with GM Tau for CFHT-Tau 4 we cannot achieve 
a fit comparable to the other objects in our sample using 
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Figure 10. SEDs and monochromatic fluxes for BDD164 where symbols have the same meanings as Figure [2] 
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Figure 11. SEDs and monochromatic fluxes for BDD304 where symbols have the same meanings as Figure [2] 
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Figure 12. SEDs and monochromatic fluxes for GM Tau where symbols have the same meanings as Figure [2] 
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Figure 13. SEDs and monochromatic fluxes for CFHT-Tau 4 where symbols have the same meanings as Figure[2] 
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any of the accreting models (see Figure 131. However, our 
derived accretion rate, logM=— 9I3J0 is consistent witli the 
literature value of —11.3 Muzerolle et al. |i2005), whilst our 



best fitting model accretion rate is significantly elevated 
(logAf=— 9). We again derive a best fitting disc mass of 
Mdisc=0.01M,. 

Whilst many stars have evidence of lower accretion rates 
as presented in Table [TJ BDD164 was found to have a lower 
accret ion rate, log M=— 11.03, and subsequently Bouy et al. 
(20081 derived a much higher accretion rate of log M=— 9.7 



using SED fitting. Therefore, within our uncertainties we are 
unable to rule out the higher accretion rates based on the 
values presented in Table [7j and more work must be done to 
ascertain whether SED fitting in general overestimates the 
accretion rates of BDD targets, or whether other methods 
systematically underestimate the accretion rates. 



5 ADOPTED PARAMETERS 

Adopting realistic uncertainties for the distance and extinc- 
tion, as oppose to fixing them, does not result in significant 
variation in the best fitting parameters (as shown in Section 
|4|. As there is no physical reason to assume all objects in 
Taurus are at 140 pc, and that the extinctions derived in 
the literature are infinitely precise, we have adopted the fits 
where these parameters are varied over a reasonable range. 
For the naked objects we have found that for most cases, 
excepting ITG 2 and CFHT-Tau 5, the SEDs are well rep- 
resented and the best fitting parameters consistent with lit- 
erature values, when fitting using the semi-empirical grid, 
(i.e. assuming negligible accretion and allowing the "model" 
or "isochronal" age to vary). Clearly, this variation in age 
simply reflects a variation in luminosity and its underlying 
causes are beyond the scope of this paper but are discussed 



in the literature (e.g. Mayne & Naylor 2008 1. For the cases of 
ITG 2 and CFHT-Tau 5 it is apparent that these objects are 
at luminosities elevated, for their masses, above the 1 Myr 
isochrone and therefore could not be fitted satisfactorily (we 
are running a grid to address this, see Section l6|. 

All of the targets with discs, BDDs, were best fit us- 
ing the accreting models of the ad-hoc grid, with the ex- 
ception of CFHT-Tau 12. This suggests that variations in 
the accretion rate and occultation by the disc are sufficient 
to account for the luminosity spreads in this population. 
Although again some objects are more luminous than the 
1 Myr isochrone namely, CM Tau and CFHT-Tau 4. 

The resulting adopted parameters for each BD and 
BDD object are included in Tables [9l [5] and |10| for naked 
and disc objects respectively. 



6 CONCLUSIONS 

We have collected a set of photometric observations for a 
sample of Taurus low-mass star and brown dwarf candi- 
dates. We then derived best fitting parameters and uncer- 
tainties for the stellar and disc properties (if required). We 
have shown that most naked brown dwarfs in Taurus can 
be fit using negligibly accreting models over a range of ages, 
with the exception of ITG 2 and CFHT-Tau 5 where the 
stars were found to be over-luminous when compared to the 



1 Myr isochrone. For targets with a disc we found satisfac- 
tory fits and derived a new list of accretion rates for all tar- 
gets using accreting models with an age of 1 Myr. The only 
exception is CFHT-Tau 12 where a different age, amongst 
small adjustments to other parameters, yielded a better fit. 
Additionally, we have shown that changes in the disc struc- 
ture can combine in complex ways to reproduce the SEDs 
of the brown dwarf and disc objects but with significantly 
different geometrical, e.g. inner disc position, properties to 
those found in the literature. As we have allowed the disc 
structure (in the ad-hoc grid) to adjust to an equilibrium 
state using a more realistic, vertical structure we favour our 
best fitting models. We have included physical mechanisms 
which are largely accepted for CTTS systems, in the ad- 
hoc grid. If the BD regime is a mere extension of CTTS 
physical regime to lower masses these mechanism will gov- 
ern the structure of BDD systems. So, by comparing these 
models with observations we are discovering limitations in 
our theoretical models, which directs improvements. On the 
other-hand fitting such systems using models with parame- 
terisations, although it leads to acceptable fits and predicted 
parameters may, obscure problems with our understanding 
of this regime. Therefore, we would encourage fitting of ob- 
servations using models derived using robust physical ingre- 
dients, as opposed to fitting to large grids constructed by 
adopting significant parameterisation, given the degenerate 
nature of SED fitting for BD and BDD objects. 

The models presented in this paper have been incorpo- 
rated into an online fitting tool. This tool currently allows 
users to select from two models grids to which observations 
can be fitted. However, the fitting tool itself is independent 
of these grids and simply requires a structured list of SEDs 
in order for observations to be fitted to any new model grid. 
We are currently running a further model grid to expand 
the parameter space covered and address some of the issues 
with the over- or under-luminous sources (when compared 
to the isochrones). It is our hope that continued develop- 
ment of the models of LMS and BD stars with discs will 
be informed by feedback from fitting tools such as the one 
presented in this paper. 

We have attempted to publish a complete description of 
our techniques on http://bd-server.astro.ex.ac.uk/ al- 
lowing users and researchers to recreate, check and compare 
to what we have done. 
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26 N.J. Mayne et al. 

APPENDIX A: DATA PROBLEMS 



We attempted initial fitting of tfie observations using tlie as- 
sociated uncertainties quoted in tlie source publications for 
the photometry. During this fitting, the optical photometry 
clearly contributed a disproportionate statistical weight to 
the fitting process. This resulted in best fitting SEDs which 
were systematically shifted away from the IR and NIR ob- 
servations. Previous studies have dealt with this issue by 
simply neglecting the shortwave, optical photometry dur- 
ing the fitting process ( Guieu et al.|[2007| . However, given 
the powerful constraint such observations can provide on 
the stellar photosphere we feel such observations must be 
retained, where possible. 



This initial fitting essentially adopted uncertainties 
based on a frequentist statistical approach, where we as- 
sume that the observation recorded is a value drawn from 
an underlying distribution. The shape of this distribution is 
determined by the properties of the measuring instrument 
and reduction or processing properties. Quoted uncertainties 
on observational data points are then the values encompass- 
ing 68% percent of the integrated (about the observed value) 
assumed or modelled Gaussian uncertainty curve. These un- 
certainties are by nature underestimates (encompassing 68% 
of the probability curve) and are to be understood as the 
probability that the observation, under the same exact con- 
ditions (including time), would be repeated. We adopt a 
Bayesian perspective whereby the uncertainties are a prop- 
erty of the model not the data. Such uncertainties must in- 
clude a component accounting for poor or incomplete model 
formulation, and a component accounting for uncertainties 
in conversion of the model to the observational domain. In 
effect we are asking which model, from our set, provides the 
best representation of the data, as oppose to, from which 
model distribution are the data most likely drawn. In prac- 
tice, it was clear from the initial fits that the quoted ob- 
servational uncertainties were insufficient to account for all 
sources of discrepancy between the model and observations, 
particularly for the optical regime. The resulting best fitting 
models were systematically shifted towards the shortwave 
length observations leading to a shift away from the longer 
wavelength points. Therefore, by adjusting the uncertainties 
iteratively with the fitting process, for the entire sample, we 
constructed a representative model of the total uncertainty 
for each photometric regime, short, mid and long wavelength 
(these are presented in Table [2|. Practically, this resulted in 
a significant widening of the quoted photometric uncertain- 
ties, which was most extreme for the optical photometry. 

In this study we have combined many disparate sources 
of photometry using different optical equipment (telescopes, 
cameras and filters) spread over several epochs, whilst the 
models are time independent. We have also combined sev- 
eral model formulations, for instance models of the stellar 
interiors, atmospheres and models of radiative transfer in 
discs. Therefore, our combined uncertainties, applied to our 
models, will have several contributions. In this section we ex- 
plore some of the underlying contributions to observational 
uncertainties in each photometric regime. 



A0.4 IRAC Photometry 

For the IRAC photometry an uncertainty of 0.2 mags was 
found to best represent the required uncertainties. Is this re- 
fiected by investigation of the data? Tables [XT] and [X2] show 
IRAC photometry for each target object, as published in 
several studies, for objects in our sample both without and 
with discs respectively. In almost all cases the magnitudes 
of the IRAC photometry vary over at most four observa- 
tions, much more than the quoted uncertainties. In some 
cases the differences could be attributed to improved re- 



ductions (e.g. Luhman et al. 2006 20101, whilst in others 
new observations have been performed (e.g. Luhman et al. 
^2006j ^Monin et al.||2010 l. In either case it is clear that the 
natural variation in such photometry is around a factor ten 
larger than the quoted uncertainties. Some of this variation 
will be caused by temporal changes in the structure of the 
associated disc, as demonstrated by the fact that the dis- 
persion of the photometry is generally less for naked stars 
(Table [Al| than those with discs ( Table [A2| |. However, there 
is still variation in naked stars in excess to the quoted un- 
certainties, which could be due to variations in the stellar 
emission. Even without natural variation in the photometric 
magnitudes it is clear that when combining several sets of 
photometric observations (even using the same instrument) 
one may expect a variation in the magnitudes larger than 
the quoted uncertainties. Indeed, a value of 0.2 for the ap- 
proximate uncertainty across the sample is not inconsistent 
with the spread of magnitudes shown in Table |A1| and |A2[ 



A0.5 JHK Photometry 

For the JHK photometry we found an uncertainty of 0.2 
mags to best represent the combined uncertainties of the 
models. However, in almost all the sources quoted in the 
caption of Table[2]the JHK magnitudes were consistent, this 
is as they are all taken from the 2MASS catalogue. However, 
one would expect temporal variation in the JHK photometry 
of ~0.2 mag over the timescale of a few days, as shown 



by Carpenter et al. (20011. Therefore, when combining the 



JHK observations with non-contemporaneous observations, 
reflecting the star-and disc properties at a different time, 
this variation must be incorporated into the uncertainties. 



A0.6 Optical Photometry 

The final regime where we must model our uncertainties is 
in the optical. Photometry in bands shortward of R and 
/ proved impossible to include in the fitting. The associ- 
ated uncertainties required to create a sensible weighting be- 
tween such observations and longer wavelength magnitudes, 
resulted in these observations providing no real constraint to 
the fit. For observations in the R and I bands uncertainties, 
again for the entire sample, of 0.8 and 0.3 mags respectively 
were found to result in the most appropriately weighted fits. 
Observations in this regime are the most uncertain due 
to the relatively cool effective temperatures of BDs and LMS 
(compared to MS stars). These objects, therefore, appear 
faint at these wavelengths observations in optical filters are 
difficult to obtain with any precision. Additionally, cool at- 
mospheres contain many more molecular species and current 
atmospheric models are known to include incomplete lists 
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Name 


[3.6] 


[4.5] [5.8] 
Magnitudes 


[8.0] 


Sources 


KPNO-Tau 4 


12.49±0.04 

12.52 

12.57 
12.56±0.02 

0.08 


12.34±0.06 

12.34 

12.37 
12.35±0.02 

0.03 


12.20±0.06 

12.14 

12.21 
12.23±0.04 

0.09 


12.11±0.06 

12.10 

12.08 
12.11±0.04 

0.03 


Luhman et al.||2006 
Cuicu ct al. (2d07| 
Monin et al.|||2010l 






Luhman et al. 


\\Mo\ 


~ A 






CFHT-Tau 15 

~ A 


13.17±0.03 

13.19 

13.24 
13.20±0.02 

0.07 


13.06±0.07 

13.16 

13.15 
13.05±0.02 

0.11 


12.95±0.06 

13.04 

13.25 
12.98±0.04 

0.30 


12.86±0.10 

13.04 

13.06 
13.00±0.07 

0.20 


Luhman et al.'('2006| 
Guicu ct al. (2007; 
Monin ct al. (20ldl 
Luhman et al.|(|20lo| 


KPNO-Tau 5 

~ A 


11.00±0.06 

11.03 

11.05 
11.05±0.02 

0.05 


10.90±0.05 

10.99 

11.02 
10.92±0.02 

0.10 


10.84±0.04 

10.90 

10.94 
10.81±0.03 

0.13 


10.79±0.04 

10.84 

10.83 
10.85±0.03 

0.06 


Luhman et al. 720061 
Guicu ct al. (2007) 
Monin ct al. (2010l 
Luhman et al.^(^2010| 


CFHT-Tau 16 

~ A 


13.15±0.06 

13.26 

12.23 
13.21±0.02 

0.11 


13.10±0.06 

13.16 

13.16 
13.11±0.02 

0.06 


12.95±0.09 

13.01 

13.04 
13.03±0.05 

0.09 


12.83±0.09 

12.97 

12.99 
13.02±0.06 

0.19 


Luhman et al.'('2006[l 
Guicu ct al. (2007l~^ 
Monin et al. (2010'i 
Luhman et al.|l,20lo| 




12.79±0.05 

12.74 

12.90 
12.85±0.02 

0.16 


12.70±0.04 
12.73 
12.75 

0.05 


12.58±0.06 

12.84 

12.72 
12.73±0.04 

0.26 


12.66±0.05 
12.60 
12.67 

0.07 


Luhman et al. 


(|2006| 


CFHT-Tau 13 


Guieu et al. (|2007| 






Monin et al. (|2010|l 


~ A 








9.92±0.04 

9.83 

9.98 
9.91±0.02 

0.14 


9.77±0.05 

9.91 

9.87 
9.79±0.02 

0.06 


9.70±0.04 

9.70 

9.76 
9.74±0.03 

0.03 


9.69±0.03 

9.69 

9.72 
9.72±0.03 


Luhman et al. 


1 2006|l 


CFHT-Tau 7 
~ A&0.15 


Guicu ctal.'(^M7l 
Monin et al. (201o| 
Luhman et al. 1(12010 






10.37±0.04 

10.48 

10.46 
10.42±0.02 

0.11 


10.21±0.05 

10.25 

10.27 
10.19±0.02 

0.08 


10.06±0.04 

10.03 

10.10 
10.08±0.03 

0.07 


10.02±0.03 

10.01 

10.07 
10.05±0.03 

0.06 


Luhman et al. 


||2006l 


CFHT-Tau 5 


Guicu ct al. (20071 
Monin et al. (201o| 






Luhman et al. 


2010 




~ A 








11.12±0.07 

11.13 

11.19 
11.19±0.02 

0.07 


11.03±0.06 

11.09 

11.12 
11.06±0.02 

0.09 


10.98±0.04 

10.96 

11.04 
11.02±0.03 

0.08 


10.94±0.04 

10.96 

10.99 
10.98±0.03 

0.05 


Luhman et al. 


(12006 




CFHT-Tau 11 


Guieu et al.|(|I 


(]07| 






Monin et al. |^ 


2010 






Luhman et al. 


^oJ 







~ A 






KPNO-Tau 9 


13.51 

13.63 

13.61±0.02 

0.12 


13.52 

13.53±0.02 

0.01 


13.67 

13.44±0.05 

0.23 


13.41 

13.48±0.08 

0.07 


Guieu et al.l (120071 
Monin et al. (201(|_^ 
Luhman et al.|(|2(J10 




~ A 




CFHT-Tau 2 


11.62±0.05 

11.54 

11.63 
11.57±0.02 

0.09 


11.38±0.05 

11.41 

11.43 
11.38±0.02 

0.05 


11.35±0.04 

11.32 

11.34 
11.37±0.03 

0.05 


11.29±0.03 

11.33 

11.33 
11.36±0.03 

0.07 


Luhman et al. ( 2006 
Guicu ct al. (2007) 
Monin et al. (2010j 
Luhman et al.|(|2010 




~ A 






11.71±0.05 

11.78 

11.79 
11.77±0.02 

0.08 


11.62±0.05 

11.69 

11.69 
11.62±0.02 

0.07 


11.54±0.04 

11.61 

11.60 
11.56±0.03 

0.07 


11.56±0.04 

11.59 

11.57 
11.54±0.03 

0.05 


Luhman et al. 


(12006 




CFHT-Tau 3 


Guieu et al.|(|I 


(J07 






Monin et al. | 


2010 






Luhman et al. 


IJ^O 







~ A 








9.53±0.04 
9.56 
9.60 

9.57±0.02 
0.07 


9.41±0.04 

9.45 

9.47 
9.42±0.02 

0.06 


9.31±0.03 

9.31 

9.37 
9.33±0.03 

0.06 


9.29±0.03 

9.28 

9.31 
9.34±0.03 

0.05 


Luhman et al.|l|20C 


§ 


ITG 2 


Guicu ct al. (2007| 
Monin et al.^(^2010 


L 






Luhman et al. 


2010 




~ A 
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Name 


[3.6] 


[4.5] 


[5.8] 


[8.0] 


Sources 






Magnitudes 








11.13±0.04 
11.14 
11.16 


10.85±0.05 
10.86 
10.88 


10.49±0.03 
10.45 
10.51 


9.80±0.03 
9.80 
9.83 


Luhman et al. 


2006| 


CFHT-Tau 9 


Guieu et al. (|2 


m\ 






Monin et al 'fSOlCi 


> 




11.11±0.02 
0.05 


10.80±0.02 
0.08 


10.48±0.03 
0.06 


9.85±0.03 
0.05 


Luhman et al. 


(2fll0| 


~ A 








13.07±0.06 


12.74±0.06 


12.36±0.08 


11.59±0.07 


Luhman et al.'('2006ll 


KPNO-Tau 6 


13.08 

13.12 

13.08±0.02 

0.05 


12.79 

12.77 

12.77±0.02 

0.05 


12.47 

12.42 

12.41±0.04 

0.11 


11.68 

11.58 

11.82±0.04 

0.23 


Guicu ct al. (2007) 
Monin et al.|l|2010| 






Luhman et al. 


\Mti\ 


~ A 








12.54±0.04 


12.23±0.03 


11.89±0.05 


11.24±0.02 


Luhman et al. 720061 
Guicu ct al. (2007) 
Monin ct al. (2010) 


KPNO-Tau 7 


12.59 
12.62 


12.27 
12.29 


11.90 
11.99 


11.22 
11.25 




12.60±0.02 


12.27±0.02 


11.93±0.03 


11.26±0.03 


Luhman et al.^(^2010| 


~ A 


0.08 


0.06 


0.10 


0.04 


CFHT-Tau 12 


10.86 

10.77±0.02 

0.09 


10.63 

10.54±0.02 

0.07 


10.34 

10.31±0.03 

0.03 


9.95 

9.93±0.03 

0.02 


Monin et al. J2010|) 




Luhman et al. 


\Mti\ 


~ A 








10.71±0.03 


10.12±0.05 


9.55±0.03 


8.86±0.03 


Luhman et al.V2006|l 


BDD399 


10.77 


10.19 


9.66 


8.91 


Guicu ct al.' (20071 


10.80 


10.21 


9.64 


8.92 


Monin ct al. (2010;) 




10.79±0.02 


10.14±0.02 


9.61±0.03 


8.91±0.03 


Luhman ct al.|l|20lo| 


~ A 


0.09 


0.09 


0.11 


0.06 




9.16±0.05 
9.25 


8.70±0.07 
8.76 


8.38±0.03 

8.38 


7.79±0.02 
7.80 


Luhman et al. 


l|2006| 


GMTau 


Guleuet al.l(<^(507|l 




9.27 


8.77 


8.43 


7.81 


Monin ct al. (2010) 




9.44±0.02 


8.95±0.02 


8.64±0.03 


7.97±0.03 


Luhman et al.|lj20lo| 


~ A 


0.28 


0.25 


0.26 


0.18 


CFHT-Tau 6 


10.66±0.05 

10.78±0.02 

0.12 


10.37±0.06 

10.37±0.02 

0.00 


9.93±0.03 

10.03±0.03 

0.10 


9.10±0.03 

9.16±0.03 

0.06 


Luhman et al. 


2006 


Luhman et al. 


2010 


~ A 








9.38±0.05 
9.48 
9.54 


8.97±0.06 
9.06 
9.07 


8.54±0.02 
8.58 
8.60 


7.79±0.03 
7.79 
7.79 


Luhman et al. 


2006 




CFHT-Tau 4 


Guicu ct al. [2 
Monin ct al. ( '^ 


(07% 
MIOJ 






9.51±0.02 


9.08±0.02 


8.62±0.03 


7.85±0.03 


Luhman et al'.H^Olo| 


~ A 


0.16 


0.11 


0.08 


0.06 




10.85±0.06 
10.83 


10.29±0.03 
10.31 


9.86±0.04 
9.86 


9.19±0.02 
9.18 


Luhman et al.|(|20 


06JI 


CFHT-Tau 8 


Guicu ct al. ( 2007 






10.23±0.02 




9.14±0.03 


Luhman ct al.^(^2oI()| 


~ A 


0.02 


0.08 


0.10 


0.07 




11.38±0.05 
11.37 
11.43 


10.85±0.04 
10.88 
10.93 


10.40±0.04 
10.43 
10.50 


9.52±0.03 
9.53 
9.54 


Luhman et al. 


2006 




BDD304 


Guieu etal.|(|2(50Tl 
Monin ct al. (2010'*) 








10.87±0.02 




9.59±0.03 


Luhman ct al.|l,201o| 


~ A 


0.06 


0.08 


0.10 


0.07 




9.48±0.05 
9.51 
9.56 


8.92±0.06 
8.99 
9.00 


8.28±0.03 
8.33 
8.36 


7.40±0.03 
7.40 
7.43 


Luhman et al. 2006 




BDD164 


Guieu et sAS^'lOOt 
Monin ct al. (2010 








9.75±0.02 


9.34±0.02 


8.87±0.03 


7.89±0.03 


Luhman et al.||i20 


0) 


~ A 


0.27 


0.32 


0.59 


0.49 





Table A2. IRAC photometry for our target objects with discs as reported by several studies. A here is the approximate range of the 
quoted magnitudes. 
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of the spectral lines, on which the opacity depends IStauffer] 
et al. (2007). However, the most significant problem with 



optical photometry is uncertainty about the natural system 
of the observations. Small changes in the filter responses, es- 
pecially for asymmetric filters such as the R filter, will lead 
to large changes in the derived fluxes, as the SED gradient 
is steep across this wavelength range. 

The R and / photometry presented in Table [2] comes 
from four sources. 



Martin et al. 



(|2001| used the CFHT f 2K 
camera and state that they observe using the Rc,Ic fil- 
ters (although the CFHT f 2K camera only has Mould R, I 
filters). They go on to say that the photometry will be ex- 
plained in further detail in Dougados et al (in preparation) 
which we assume refers to Dougados et al. (200f I, although 



no further explanation is provided. Therefore, it is unclear 
whether these photometric observations are in the CFHT 
12K-Mould R, I system or have, in some way, been trans- 
formed to the Cousins system. The transformation between 
photometric systems is usually done using observation of 



Landolt standards ( Landolt f992) . However, for BDs and 
LMS such transforms are likely to be extremely unreliable 
as the SED has a steep gradient across the filter, meaning 
uncertainties in the filter shape will lead to large inaccura- 
cies in the photometry. Briceiio et al. ( 2002 1 observed us- 



ing the KPNO / filters and calibrated this to the Cousins 



system using Landolt (19921 standards. Luhman (20041 ob- 



serves using an I filter on the four-shooter camera at the 
Fred Lawrence Whipple Observatory, and employs conver- 
sions using " 



Guieu et al.] ( |2007 K 
CFHT 12K R and / 



Landolt ( 1992 I standards. The final source is the 



series of papers Guieu et al. (20051, Guieu et al. (20061 and 



Guieu et al. 



( 2005 I observe using the 
and the Mega Prime r' and i' filters, 
and the calibration is not explained, suggesting that these 
measurements are presented in the natural systems in which 
they are taken. These data are extended by observations pre- 



sented in Guieu et al. ( 2006 1 where observations, again using 
the CFHT 12K and Mega Prime systems, were combined us- 
ing overlapping fields, resulting in a photometric magnitudes 
in the Mega Prime r' and i' system, however the magnitudes 
are not explicitly presented. Finally, R and I magnitudes 
Guieu et al. (20071 where they state: "All R and 



appear m 

I data have been transformed to the CFHT12k camera's 



Cousins system (see Guieu et al.||2006 Briceno et al. 2002| ) 
before conversion to absolute fiuxes." However, the magni- 



tudes presented in Guieu et al. ( 2005 1 match those presented 
(for common objects) in Guieu et al. (20071. Therefore, as 



no transformation or calibration to Cousins is explained in 

([2006) 



either Guieu et al. ( 2005 1 or Guieu et al 



we assume 



(20071 are 



that the magnitudes presented in Guieu et al. 
actually in the Mega Prime r' , i' system. 

Table [A3] shows sources and photometry for the targets 
where we have duplicate I band measurements (all stated to 
be in the Cousins I photometric system). It is clear that for 
stars both with and without discs large variation, as much as 
0.23 mags, is apparent in the / band photometry. It is rea- 
sonable to assume that the predictive power of our models 
will decrease towards shorter wavelengths as the dust opac- 
ity becomes more important and the slope of the SED across 
the filter responses leads to an exaggeration of slight errors 
in the adopted response. Therefore, we conclude that a min- 
imum of 0.2 mags could be contributed to the uncertainty 
budget due to variation in the photometry. 



Name 



R 



Source 



Sources without infrared excess 



KPNO-Tau 4 
KPNO-Tau 5 
CFHT-Tau 7 
CFHT-Tau 5 
CFHT-Tau 2 
CFHT-Tau 3 
ITG 2 



20.54 
19.10 
16.63 
23.37 
20.21 
20.33 
20.21 



Guieu et al.|j2006| l, Guieu et al. (|2007l 
Guieu et al 



Guieu ct al (2005), Guieu et al 
Guieu ct al. ( 2 005 J 
Martin ct al. ( 2Ml I 
Martin et al.iiOTTi 



I 2007 1 



Guieu et al. 



il|20PI 



Sources without infrared excess 



KPNO-Tau 6 
BDD399 
CFHT-Tau 6 
CFHT-Tau 4 
CFHT-Tau 8 



20.56 
20.33 
18.39 
19.10 
19.27 



Guieu et al. 



Guieu et al. 



120051, 



mm 



Guieu et al. 



I 2007 1 



Guieu ct al. (2005) 
Martin et al.|||2001|) 



Guieu et al. 



ilj2001|) 



Guieu et al. I 2007 1 



Table A4. A table showing the -R band photometry and sources. 



Table |A4| lists the sources of the R band photometry. 
For the R band photometry there are no duplicates (where 
the photometry is from different underlying sources), so we 
cannot assess the uncertainty of the transforms applied (as 
all of the photometry is stated as having been transformed 
to the Cousins system). However, we can assume that the 
uncertainty will be significantly larger than that found for 
the / band photometry due to the shapes of the R band 
filters. Figure Al shows the filter responses (normalised to 



the same scale) for the Mega Prime, CFHT and Cousins R 
and / filters. The shapes of the I band filters are roughly 
consistent with the assumption of a fiat response. For the 
R band however there exists a longer tail (extending to 0.0 
^m) in the Cousins filter ( |Bessell|2005[ ). This will exacerbate 
problems with transformations between these filters as the 
slope of the SED is a sharp function of wavelength in this 
region (see Figure [AT| . 

Transformations are generally performed using |Lan-| 
dolt ( 1992 ) standards, which are unreddened main-sequence 



(MS) stars in a limited colour range. To demonstrate why 
such transforms are unreliable for BD stars we have per- 
formed a simple experiment to construct a simulated ver- 
sion of a transformation between the relevant photometric 
systems. In practice transformations are made using a con- 
version equation as a function of colour. Table [AS] shows the 
simulated R and / band photometry in the Cousins, CFHT 
12K and Mega Prime systems for three unreddenedj MS 
stars at colours _Rc — Ic ~0, 1.2 and a pre-MS BD at 2.3 
behind 2 magnitudes of extinction (typical for the targets 
in Taurus). The simulated photometry was created using 
the models of|Siess et al.l(|2000[),|Chabrier et al.| ( ;2000J and 



Hauschildt et al. ( 1999). Briceno et al. (2002), for example. 



used the standard fields SA 95 and SA 101 to calibrate their 
photometry to the Cousins system. These standard fields 
contain MS stars covering a range of colours from ~0.0 to 
1.2. Whereas, the pre-MS stars of colour 2.3 is typical of 
the bluest BD stars in Taurus. 

Table IA5I shows that if we construct transformation 



^ at 140 pc, although distance is unimportant this has been set 
to the distance of Taurus for illustrative purposes 
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Name 



I Sources 

Magnitudes 



CFHT-Tau 2 



KPNO-Tau 6 

~ A 



CFHT-Tau 6 

~ A 



CFHT-Tau 4 

~ A 



CFHT-Tau ! 

~ A 



Sources without infrared excess 



16.81 
16.69 
0.12 



Briceno et al. I 2002 I 
Martin et al. |200l| 



Guieu et al. 



I 2005 I 



Guieu et al. 



I 2007 1 





16.88 
16.77 
0.11 


Bricefio et al. 


jjnn.'[,. 


Guieu et al. 


2007|l 


CFHT-Tau 3 


JVLartin et al. 


m\\ 






~ A 









Sources with infrared excess 



17.90 
17.80 
0.10 



Guieu et al. 



Briceiio et al 



( 2005 1^ Guieu et al. 
1^2 



I 2007 1 



15.40 
15.63 
0.23 



Guieu et al. 



Luhn 



(20051 



^Oli 



15.78 
15.64 
0.14 



Luhman 



Martin et 



1 2004 1 

|2001 l 



il 



16.43 
16.24 
0.19 



Guieu et al. 



Luhn 



(20051 



"^^^ 



Guieu et al. 



I 2007 1 



Table A3. Duplicate photometry for which differences are above 0.05 mag 



Rc-Ic 


Band 


Cousins 


CFHT 12K 


Mega Prime 






pre-MS, 


Av =2 




~2.3 


R 

I 

R-I 


19.63 
17.31 
2.32 


20.31 
17.08 
3.23 


20.87 
17.71 
3.16 






MS, A 


v=0 




~0 


R 

I 

R-I 


6.35 
6.37 
-0.02 


6.36 
6.38 
-0.02 


6.37 
6.37 
0.00 


~ 1.2 


R 

I 

R-I 


10.93 
9.73 
1.20 


11.08 
9.68 
1.40 


11.26 
9.95 
1.31 



Table A5. Table showing the simulated photometry for MS stars 
across a range of colours and a pre-MS BD typical of the Taurus 
targets. 



equations (assuming a linear relationship in colour between 
the instrumental and calibrated magnitude) using the MS 
stars in the colour range usually found in the Landolt ( 1992 1 



standard fields we derive (assuming a negligible airmass 
component) 

Ic = /cFHT - 0.009 + 0.042 X (i?cFHT - /cfht) (Al) 

Ic = /Mega + 0.000 - 0. 168 X (i?Mega - hicg.) ( A2) 

Re = i?CFHT - 0.012 - 0.099 X (i?cFHT - /cfht) (A3) 

Re = /?Mcga - 0.020 - 0.252 X (/^Moga - /Mega). (A4) 

These relationships would then be used to extrapolate out to 
the required colour for transformation of Taurus BDs. The 



Band 


CFHT 12K 


Mega Prime 


-'transforms 
-'simulated 


+0.127 
+0.230 


-0.531 
-0.400 


^transforms 

-'^simulated 


-0.332 
-0.680 


-0.816 
-1.240 



Table A6. Table showing the A magnitudes predicted from the 
transforms equations which simulate the use of [Landolt] | |1992| 
standards, and the actual A magnitudes for a typical simulated 
Taurus object. 



results of evaluating the transform equations and comparing 
to the values derived from the simulated observations are 
shown in Table lA6l 



Table A6 shows that if Landolt ( 1992 1 transforms are 



used we would expect the transforms to be incorrect by 
around ~0.1 for the / band and ~0.4-0.5 for the R band 
for pre-MS BDs. This is caused by the shape of the SED 
changing across the filter responses as one moves from MS 
stars to BDs. The change in SED between BDs and MS is 
driven by a change in the surface gravity of these objects, 
~5.25 for a MS star and 3.79 for a pre-MS star at similar 
colours. Figure \Kl\ shows the normalised filter responses for 
the / and R bands of the Cousins, CFHT 12K and Mega 
Prime filter systems, overlaid with SEDs of a MS BD and 
pre-MS BD (arbitrarily scaled) for illustrative purposes. 

Figure [XT] shows that the area under the filter response 
is similar for the / band filters but the R band Cousins filters 
shows a longer wavelength tail. The overlaid SEDs then show 
that for the MS (greater flux, or higher line) and the pre-MS 
the long wavelength tail of the Cousins filter becomes signif- 
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Figure Al. The normalised filters responses as a function of 
wavelength for the / (left panel) and R (right panel) filters 
of the Johnson-Cousins system (green solid line |Bessell||2005l l, 
Mega Prime (blue solid lines, r', i' of http://cfht .hawaii.edu/ 
Instruments/Filters/) and CFHT 12K (red solid line, Mould R 
and /, http://cfht.hawaii.edu/Instrimients/Filters/). Over- 
laid, normalised to arbitrary units, are the SEDs (dashed lines) 
of a MS (unreddenned) and pre-MS {Ay =0) star at a colour 
typical of a Taurus object {R — 7=2.3). 



leant. This leads to the problems with transforms illustrated 
by Table [XG] Therefore, for the / band photometry the com- 
bination of variations in the photometry and transformation 
between photometric systems, leads to an uncertainty bud- 
get of ~ 0.3, matching the value we have adopted for the 
sample. For the R band photometry, uncertainties caused 
by transformations, combined with variation in the photo- 
metric magnitudes, lead to uncertainties of ~0.7. However, 
during our fitting process we found uncertainties of 0.8 were 
required, suggesting that some further source of uncertainty 
may be required, for instance the missing opacity sources in 
the BD atmosphere models ( Stauffer et al.|2007 l. 

As discussed previously other authors have also tried 



to perform fitting using the NOMAD catalogue (Zacharias 
et al.|2005 l, which is derived from photographic plate mea- 
surements leading to extreme uncertainty in their conversion 
to Johnson-Cousins magnitudes, therefore as we are unsure 
of how to transform this photometry correctly we have not 
included these measurements. 

Therefore, for BD objects, which have low effective tem- 
peratures (compared to MS stars), photometric transforma- 
tions, within the optical regime, are likely to produce sig- 
nificant uncertainties in the resulting magnitudes, which are 
exacerbated due to the shape of the R filter. 



APPENDIX B: CALIBRATION AND 
SIMULATED PHOTOMETRY 

Initially, our simulated stellar photospheres were diluted by 
the factor (/?, /distance)^ to a distance of 10 pc to produce 
SEDs in absolute flux. These SEDs then represented simu- 
lations of 'naked' stars, which, as the emission is isotropic 
are inclination independent. For those systems where the 



disc has been modelled similar absolute flux SEDs were also 
constructed at the inclinations detailed in Table [S] 



Bl Photometric Systems and Calibration 

To allow comparison and fltting to real observations these 
SEDs must be adjusted to account for the distance and ex- 
tinction to a given target. We chose an extinction system 



consistent with Robitaille et al. (20071, using the same ex- 



tinction curve ( Indebetouw et al. 2005 1 interpolated onto our 



wavelength sampling. This aflowed the construction of SEDs 
for stars with and without discs at any distance and extinc- 



tion (using Fa] 



Fab. X 10 



(-OAA^) 



10 
^ distance{pc) ' 



where 



Fapp is apparent flux, Fabs absolute flux and Ax the extinc- 
tion in the required band) , for comparison to observed spec- 
tra. However, most observations of such systems are photo- 
metric and presented as either monochromatic fluxes or in a 
specifled magnitude system. Therefore, to enable fitting we 
derive these photometric quantities. 

To derive magnitudes one requires the underlying fil- 
ter responses of the instrument, the method by which the 
instrument sums flux (i.e. counting photons or energy) and 
the assumed zero points. 

The magnitudes are deflned as 



M. 



-2.51ogj^Q(A''^°"" '') — Zeropoint, 



(Bl) 



where counts is usually the number of electrons, or in the 
case of IRAC where counts=e^/Gain, as the gain is sig- 
niflcant (~3-4). Subsequent integration across the filter is 



performed as follows (as described in Bessell et al. 1998 1 



For photon counting, i.e. responses in electrons/photon. 



iV" 



f^^Rx{X)Fx{\)d\ 



(B2) 



and for energy summing, i.e. responses electrons/unit en- 
ergy, 



iV" 



/ 



Rx{\)Fx{\)dX 



(B3) 



The derivation of monochromatic fiuxes is explained in Ap- 



pendix A of Robitaille et al. (20071. The summed flux is 



still calculated, as with magnitudes, but this flux is then 
redistributed across the filter using an underlying (well 
chosen) assumption as to the SEDs shape (for instance 
10 000 K Black body for the Spttzer MIPS fiuxes). Finally, 
the "monochromatic" fiux is then quoted as the value of this 
redistributed fiux curve at the central wavelength of the fil- 
ter. After integration using equations |B3| or [B2l an additional 
normalisation factor is required. For a monochromatic flux. 



/t/(A)FA[actual]i?A(A)dA _ /[/(A)FA[assumed]i?A(A) dA 



jRx{X)dX 



jRx{X)dX 



(B4) 



where U{X) is 1 if the response is in e~/unit energy or X/{hx 
c) if the response is in e~ /photon. Then, as the assumed 
shape is known, 



XFx[assumed] = XoFxg[quoted]. 



(B5) 



Therefore, to correctly calibrate magnitude and monochro- 
matic flux derivations one require the fllter responses, as- 
sumed SED shape and photometric zero point. A summary 
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of all of the calibration information for our derived photome- 
try is available onlinajand is presented in Table |Bl[ Present- 
ing this information is vital as it means that researchers can 
recreate what we have done and compare quantities derived 
from the models with the same observational quantities re- 
gardless of the system chosen. 



B2 Fitting Photometric Measurements 

The derivation of magnitudes and monochromatic fluxes is 
performed during fitting on the apparent SEDs. At this stage 
the reader may ask the question why do we spend the extra 
computational time to derive the magnitudes and fluxes di- 
rectly (at run time) from the apparent SEDs, as oppose to 
simply adjusting the photometry from the absolute model by 
the extinction and distance. Additionally, the reader might 
question why we derive the monochromatic fluxes so rigor- 
ously, as oppose to simply interpolating for a flux point at 
the assumed wavelength on the SED. 

Previous fitting, in general, has been performed using 
these two shortcuts to provide a much faster fitting rou- 
tine. However, these sacrifices imply major assumptions. Ei- 
ther, the filters are assumed to be infinitesimally thin, or 
the SED is assumed to be a smooth curve which matches 
the assumed SED shape across the filter as well as the filter 
response curves being assumed to be flat across the band- 
pass. In most, but not all, cases the uncertainties in the 
observations and data are large enough to blur any errors 
made due to the breakdown of these assumptions. 

This breakdown is apparent for objects under heavy 
extinction but is perhaps best demonstrated by the use of 
optical photometry for BDs. In this regime the SED is a 
steep function of wavelength. 

Most studies will convert optical magnitudes to so- 
called monochromatic fluxes to provide better constraint of 
the stellar parameters (e.g. Guieu et al.|2007 Luhman et al 
2007 Bouy et al. 2008 1. This conversion is usually performed 
by calculating a summed flux from the magnitude (i.e. re- 
versing equation BIl, then assuming some underlying SED 



Apparent (D=140pc Ay=10) 



shape (usually those presented in Table Bl \ and quoting the 
flux at the filters central wavelength. These fluxes are then 
fit directly to interpolated values of the flux at the same 
wavelength. 

Figure |BI| shows as solid lines the SEDs for a typi- 
cal brown dwarf and disc (BDD, hereafter) model from our 
model grid (A model from the ad-hoc grid: A/* = O.O4M0, 
Age=lMyr, negligible accretion, a rotation period of 0.5 
days, a disc mass of O.OIM* a disc outer radius of 100 AU and 
a=2.15), at 10 pc (left panel) and 140 pc behind an extinc- 
tion of lOAv (right panel). The monochromatic fluxes de- 
rived rigorously for each SED are shown for the CFHT Mega 
Prime u*, g', r' and i' filter bands (as crosses). If one applies 
the extinction to the monochromatic fluxes directly there re- 
sults a significant difference in fiux from those derived from 
the reddened SED. Additionally, the derived monochromatic 
fiuxes do not lie directly on the SED, in either case. Whilst, 
the deviation is small given the gradient of the SED in these 
region the affect on the fitting statistic is significant. Es- 
sentially, if we do not derive monochromatic fluxes properly 



(http : //bd- server . astro . ex . ac . uk/calibratlon) 




2 0.4 0.6 0.8 1.0 



Figure Bl. SEDs (solid lines) in A (/xm), XFx 
{ergs s~^cm~^ A~^), at lOpc (left panel) and 140pc be- 
hind an extinction of Ay =10 (right panel). The crosses denote 
the derived monochromatic fluxes for each SED and the circles on 
the right panel show the fluxes resulting from simply reddening 
the absolute monochromatic fluxes. 



from the underlying SEDs, and furthermore, simply treat 
them as points along the SED we are not comparing two 
equal quantities. Therefore, for our fitting we only compare 
observed monochromatic fiuxes with those derived from the 
underlying SEDs in the same way. 

This problem is exacerbated in the optical regime as, 
due to the steep gradient across the wavelength space, the 
monochromatic fiuxes in this region are very powerful con- 
straints. In addition the photometric systems and conver- 
sions of these fluxes from magnitude space are often un- 
clear or unreliable. For instance several studies have utilised 



the photometry of Zacharias et al. (20051 or the NOMAD 



survey and converted the optical magnitudes to monochro- 
matic fluxes. However, these magnitudes are often derived 
from photographic plates and it is not obvious how to con- 
vert these correctly to comparable monochromatic fluxes. In 
fact transformations between magnitude systems are usually 
derived using observations of main-sequence stars and then 
applied, without adjustment, to LMS and BDs (this is dis- 
cussed more in Section [A|. 

A further complication for overlaying monochromatic 
fiuxes on spectra (if one does not derive the monochromatic 
fiuxes rigorously) is the quoted central wavelength. In real- 
ity the effective wavelength of a filter shifts position with 



the input spectrum as shown by Bessell et al. ( 1998 1 and is 
defined as, 



Ae 



^ \Fx[actual]Rx{\) d\ 
J Fx[actual]Rx{\) d\ ' 



(B6) 



For a typical spectrum of a BD the shift in wavelength is 
small for filters at longer wavelengths (typically ~0.01/im). 
However, for shorter wavelength filters, where the photo- 
spheric fiux changes rapidly across the filter, the shift can 
be as large as ~0.1/im. If one fits observed monochromatic 
fiuxes to monochromatic fluxes derived for an SED in the 
same fashion then the effective wavelength of the fllter does 
not affect the difference (it is just a direct comparison of 
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fluxes). However, if tlie observed monocliromatic flux is fit- 
ted to a simple point lying on the SED, for optical colours, 
then the shift in effective wavelength may well be signif- 
icant, again exacerbated by the steep rise in stellar flux 
across the optical filters. Essentially, this leads to flux points 
in the optical being shifted away from the expected posi- 
tion. Therefore, one must be careful when overlaying optical 
"monochromatic fluxes" derived from optical magnitudes on 
model SEDs. This is discussed and an example presented in 
Section |4l 

As the uncertainties in the models are large most pre- 
vious studies have not required this level of sophistication 
in the derivation of photometric values. However, as this ef- 
fect becomes particularly prominent in the optical the SED 
curve and monochromatic flux points are often dropped in 
this regime (see, for example,the incomplete SEDs at lower 
wavelengths of Guieu et al.|2007l. As, using this fitting tool. 



we are now able to precisely model observations in the op- 
tical we no longer have this limitation and, additionally, as 
models improve have the software to test their predictions. 
Simply stated to improve the models we must ensure we are 
comparing them, on an equal footing, to observations. 

B3 Other Tools 

In addition to the fitting tool several other tools are avail- 
able at http://bd-server.astro.ex.ac.uk/ enabling users 
to explore the model grids. A browsing and downloading 
tool is available to select, view and retrieve, SEDs and flux 
and magnitude information for any model within any grid 
(http : //bd- server . astro . ex . ac . uk/plot s) . Additionally, 
we have constructed model isochrones and mass tracks us- 
ing our grids and these can be retrieved from http : // 
bd-server .astro.ex.ac.uk/isochrones. Finally, we have 
also made available an interpolation tool allowing users to 
produce SEDs, monochromatic fluxes and magnitudes for 
models within, but not covered by, our parameter space at 
a user entered distance and extinction. 
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Filter Set 



Response SED Shape Zero Magnitude 

Units (o-/?) (Assumed) Flux, F,(Mo) (Jy) 



Filter 



Zero 
points 



Central A 



Width 



Johnson-Cousins (1) & (2) Photons 



Flat(^) Vega(2' 



U 

B 

V 

R 

I 

J 

H 

K 



-14.18 
-15.30 
-14.86 
-15.11 
-14.58 
-14.26 
-13.49 
-13.09 



0.36 

0.46 

0.585 

0.71 

0.81 

1.23 

1.62 

2.18 



0.10 
0.18 
0.21 
0.30 
0.20 
0.36 
0.44 
0.56 



Tycho (3) & (2) 



Photons 



Flat 



Vega 



Vt 
Bt 



-15.09 
-14.94 



0.527 
0.422 



0.21 
0.145 



BesseU (4) & (5) 



Photons 



Flat 



Vega 



U 

B 

V 

R 

I 

J 

H 

K 

L 



-14.23 
-15.32 
-14.86 
-15.09 
-14.57 
-13.92 
-13.65 
-12.71 
-11.74 



0.37 
0.44 
0.55 
0.64 
0.80 
1.22 
1.63 
2.19 
3.45 



0.11 
0.18 
0.21 
0.29 
0.20 
0.36 
0.96 
0.36 
0.72 



SDSS (6) & (7) 



Photons 



Flat 



i^^(Mo)=3631 



-12.43 
-14.22 
-14.22 
-13.78 
-11.84 



0.3555 

0.473 

0.6261 

0.7672 

0.9097 



UKIRT (14), (12) & (13) Photons 



Flat 



Vega 



Z 

Y 
J 
H 
K 



-13.84 
-13.46 
-14.01 
-13.29 
-12.82 



0.8775 

1.02 

1.28 

1.65 

2.2 



0.0599 
0.1379 
0.1382 
0.1535 
0.137 









F,(Mo)=1594 


J 


15.73 


1.235 


0.162 


2MASS (8), (9) & (10) 


Photons 


Flat 


F,(Mo)=1024 


H 


16.37 


1.662 


0.251 








F„(Mo)=666.7 


Ks 


17.38 


2.159 


0.262 










J 


-13.48 


1.24 


0.16 


MKO (11), (12) & (13) 


Photons 


Flat 


Vega 


H 


-13.46 


1.65 


0.29 










K 


-12.64 


2.2 


0.34 



0.095 

0.1 

0.42 

0.45 

0.66 



IRAS (15), (16) & (17) 



Energy 



Flat 



F,{Mo) 
F^Mo) 
F.{Mo) 



=28.3 
=6.73 
=1.19 
=0.43 



12 
25 
60 
100 



21.00 
23.55 
26.24 
28.41 



12.0 
25.0 
60.0 
100.0 



6.5 
11.0 
40.0 
37.0 



SCUBA^'') (18) & (19) 



Energy 



lOOOOK 
Black Body 



N/A(3) 



450WB 
850WB 



443.0 
863.0 



138.92 
230.636 



IRAC (20) & (21) 



Photons 



Flat 



F.iMo) 
F.(Mo) 
F^Mq) 
F,{Mq) 



=280.9 
=179.7 
=115.0 
=64.9 



3.6 

4.5 
5.8 
8.0 



-10.15 
-9.70 
-8.07 
-8.55 



3.55 

4.493 

5.731 

7.872 



V 



-14.86 0.5374 



CFHT 12K (Mould) (31) Photons 



Flat 



Vega 



1.2 
1.56 
2.27 
4.69 







lOOOOK 


F.(Mo)= 


=7.17 


24 


25.97 


23.68 


14.20 


MIPS (22) 


Energy 


Black 


body 


F.(Mo)= 


=0.778 


70 


29.62 


71.42 


61.06 










F,{Mo)-- 


=0.160 


160 


32.48 


155.9 


99.83 














J 


-13.73 


1.25 


0.24 


CIT (23) & (24) 


Photons 


Flat 




Vega 




H 


-13.44 
-12.97 


1.65 
2.215 


0.30 
0.41 














blue 


-6.00 


70 


90.5 


PACS (25), (26) & (27) 


Photons 


Flat 




Vega 




green 
red 


-5.26 
-4.27 


100 
160 


109 
391.25 














250 


-3.16 


250 


176.30 


SPIRE (28), (29) & (30) 


Photons 


Flat 




Vega 




350 
500 


1.87 

o.o(^' 


350 

500 


281.26 
290.24 



0.0974 
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R 


-14.69 


0.6581 


0.1251 


I 


-14.78 


0.8223 


0.2164 


u* 


-14.33 


0.3743 


0.0758 


g' 


-15.49 


0.4872 


0.1455 


r' 


-14.76 


0.6282 


0.1219 


i' 


-14.61 


0.7776 


0.1508 


z' 


-15.07 


1.1702 


0.6868 



CFHT Mega Prime (31) Photons 



Flat 



Vega 



Table Bl: A table describing the calibration information for each photometric band used. Sources are given as numbers in brack- 
ets: (l)|Johnson| ( |1966[ ), (2) |Bessell|(2005[ ), (3) |BesseIl| (|2000') , (4)'Besse ll fc Brett|(|1988[), (5 ) |Bessell et al.|(|1998[), (6)|Fukugita 
|et al.| ( |1996[ ), (7) http : //www . sdss . org/dr 1/algor ithms/f luxcal . html, (8)|Cohen et al.](|2003[), (9) [Skrutskie et al.|(|2006f7 
(10) http://www.ipac.caltech.edu/2mass/releases/allsky/doc/sec6_4a.html, (11) Tokunaga et al. (20021, (12) Simons 



& Tokunaga (20021, (13) http://www.jach.hawaii.edU/UKIRT/instruments/wfcam/user_guide/description.html#Tabl. 1, 
(14) |Hawarden et aL] ( |2001[ ), (15) [Neugebauer et al.| ( |1984'1 ), (16) http://irsa.ipac.caltech.edu/IRASdocs/archives/ 
spectral_resp.html, (17) http://irsa.ipac.caltech.edu/IRASdocs/exp.sup/ch6/C2a .html, (18) [Holland et al. (19991, 
(19) http://www.jach.hawaii.edu/JCMT/continuuiii/background/background.html, (20) Reach et al. (20051, (21) http:// 
ir sa . ipac . caltech . edu/data/SPITZER/docs/irac/, (22) http : //irsa . ipac . caltech . edu/data/SPITZER/docs/mips/, (23) 
Elias et al.] (|1982|), (24) [Stephens fc Leggett] (|2004[), (25) [Poglitsch et al.| (|2008[), (26) [Poglitsch et al.| (|2010[), (27) 



http : //herschel . esac . esa . int/Docs/PACS/html/pacs_om . html, (28) |Griffin et al.| ( |2008[ ), (29) [Griffin et al.| ( p)10| |, 
(30) http: //herschel. esac. esa. int/Docs/SPIRE/html/spire_om.html#xl-880005. 2. 3 and (31) http: //cfht .hawaii . 
edu/Instruments/Filters/. Notes are superscript numbers in bracket: (1) A Constant flux across the filter. (2) The Vega 
spectrum was a NextGen atmosphere model with roff=9950K and log((7)=3.95, scaled to the Vega flux of Hayes (19851. (4) 
SCUBA 2 filters approximately equal to SCUBA filters. (5) Wavelength longer than the maximum available in the Vega 
spectrum. 
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